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ABSTRACT 

We present results from a multiwavelength IR-to-X-ray campaign of the infrared bright 
^ ; (but highly optical-UV extincted) QSO IRAS 13349+2438 obtained with the Chandra High 

Energy Transmission Grating Spectrometer (HETGS), the Hubble Space Telescope Space 
Telescope Imaging Spectrograph (STIS), the Hobby-Eberly Telescope (HET) 8-meter, and 
the Spitzer Infrared Spectrometer (IRS). Based on HET optical spectra of [O III], we refine 
> . the redshift of IRAS 13349 to be z = 0.10853. The weakness of the [O III] in combina- 

OO ' tion with strong Fe II in the HET spectra reveal extreme Eigenvector- 1 characteristics in 

, IRAS 13349, but the 2468 km s^^ width of the H/3 line argues against a narrow-line Seyfert 1 

classification; on average, IR, optical and UV spectra show IRAS 13349 to be a typical QSO. 
Independent estimates based on the H/3 line width and fits to the IRAS 13349 SED both give 
a black hole mass of A/bh = 10^ Mq. The heavily reddened STIS UV spectra reveal for the 
' first time blue-shifted absorption from Ly a, N V and 4, with components at systemic veloci- 

I ties of —950 km s^^ and —75 km s^^ . The higher velocity UV lines are coincident with the 

lower-ionisation (^ ^ 1.6) WA-1 warm absorber lines seen in the X-rays with the HETGS. 
In addition, a ^ 3.4 WA-2 is also required by the data, while a ^ ^ 3 WA-3 is predicted 
by theory, and seen at less significance; all detected X-ray absorption lines are blueshifted 
by ~' 700 — 900kms~^ . Theoretical models comparing different ionising SEDs reveal that 
^ I including the UV (i.e., the accretion disc) as part of the ionising continuum has strong impli- 

- - cations for the conclusions one would draw about the thermodynamic stability of the warm 

absorber Specific to IRAS 13349, we find that an X-ray-UV ionising SED favors a contin- 
uous distribution of ionisation states in a smooth flow (this paper), versus discrete clouds in 
pressure equilibrium (previous work by other authors). Direct detections of dust are seen in 
both the IR and X-rays. We see weak PAH emission at 7.7 ^m and 11.3 ftm which may also 
be blended with forsterite, and 10 fim and 18 fj,m silicate emission, as well as an Fe L edge at 
700 eV indicative of iron-base dust with a dust-to-gas ratio > 90%. We develop a geometri- 
cal model in which we view the nuclear regions of the QSO along a line of sight that passes 
through the upper atmosphere of an obscuring torus. This sight line is largely transparent in X- 
rays since the gas is ionised, but it is completely obscured by dust that blocks a direct view of 
the UV/optical emission region. In the context of our model, 20% of the intrinsic UV/optical 
continuum is scattered into our sight line by the far wall of an obscuring torus. An additional 
2.4% of the direct light, which likely dominates the UV emission, is Thomson-scattered into 
our line-of-sight by another off -plane component of highly ionized gas. 

Key words: galaxies: active - galaxies — Seyfert: individual: (IRAS 13349+2439) - galax- 
ies: warm absorber - infrared, optical, ultraviolet. X-rays - observatories: Chandra, HST, 
© 0000 RAS Spitzer, HET — instruments — HETGS, STIS, HET, IRS 



2 Lee et al. 



1 INTRODUCTION 

IRAS 13349+24 38, hereafter "IRAS 13349", (z = 0.10764 - 
iKim etalj[l995l : updated here to z = 0.108530 based on the 
3x higher resolution HEX spectra reported in this paper) is a 
prototype infrared-luminous q uasar with a bolomet ric luminos- 
ity of > 2 X 10'"' erg s~^ l lBeichman et"ai] Il986l) . Images of 
the host galaxy and nearby environment show the galaxy to be 
spiral-like , with a possible companion at ~ 5" along the minor 
axis (e.g., iHutchings & McClurdll990l) . Evidence for tidal struc- 
ture suggests that the object may have interacted with the com- 
panion, and this could supply gas and dust to the nucleus, fueling 
quasar activity and enhancing nuclear obscuration and scattering. 
Indeed, IRAS 13349 has a broad-emission line optical spectrum 
that becomes heavily reddened at shorter wavelengths, and exhibits 
high optical continuum and emission-line polarisation jWills et al] 
I1992I hereafter W92). The observed polarisation rises strongly to- 
ward shorter wavelengths, but the optical polarised flux spectrum 
is indistinguishable from a typical, unreddened quasar (Hines et al. 
2001). These polarisation properties indicate that observers see the 
quasar's nucleus via both a direct, but attenua ted, light path and 
a scattered light path (W92; iHines et al]|200lh . In addition to its 
high luminosity and spectral variability, IRAS 13349 also exhibits 
strong Eigenvector- 1 characteri stics (strong optical Fe 1 1 emission 
and weak [O III] relative to Hg: lBoroson & Greenlll992h . It is typ- 
ically identified to be radio-quiet, altho ugh weak 4.87 mjy radio 
emissi on has been reported at 6 GHz bv lLaurent-Muehleisen et aU 
( I1997I) . 

The observed high X-ray flux and large-amplitude X-ray vari- 
ability indicate that a large fraction of the X-ray emission is seen 
via a direct, rather than a scattered, path. As such, IRAS 13349 
has been the subject of a number of X-ray studies that have 
helped to clarify the physical processes in the inn er regions of the 
quasar nucleus. Using XMM-Newton EPIC data, iLonginotti et al.l 
( l2003h made a detailed analysis of the IRAS 13349 ionised re- 
flection spectrum and relativistic Fe emission line. At lower en- 
ergies . X-ray studies based on the ROSAT PSPC jSrandt et al.l 
1 19961) and ASCA l^andt et al. 1997), combined with optical/near- 
infrared extinction estimates argue for obscuration by dusty, 
ionised gas. Studies with modem day high spectra l resolution in- 
struments on board XMM-Newton jSako et al]|200ll) . and Chandra 
jHolczer et al .l2007l : lBehaill2009l : also this paper) reveal additional 
complexity in the absorbers (most notable an unresolved transition 
array of 2p-3d inner-shell absorption by iron M-shell ions, dubbed 
the UTA by its discoverers, Sako et al.), and allow direct measure- 
ments of the dust composition (this paper) in the host galaxy. 

In this paper, we present a comprehensive analysis of the ab- 
sorber properties of this quasar, based on our high spectral resolu- 
tion, multi-wavelength campaign, involving X-ray (with the Chan- 
dra High Energy Transmission Grating Spectrometer; HETGS), ul- 
traviolet (with HST STIS-MAMA), optical (HET: Hobby-Eberly 
Telescope 8-m), and infrared (Spitzer IRS) observations. The 
Chandra and HST observations are simultaneous, HET near simul- 
taneous, and Spitzer IRS taken 1.25 years later as part of the GTO 
program of G. Rieke. These high-quality data allow us to address, 
in detail for a specific luminous system, several issues of broader 
importance in active-galaxy absorption studies. These include the 
apparent presence of dust in some ionized absorbers and its im- 
plications for interpreting observations; the relation of ionized ab- 
sorbers to other nuclear components, including accretion disks, tori, 
and scattering material; the relation of the spectral energy distribu- 



Table 1. Chandra Observations 



Date 


Start Time 


End Time 


Useable Time 




(GMT) 


(GMT) 


(ks) 


2004-02-22 


06:42:06 






2004-02-23 




12:15:15 


160 


2004-02-24 


13:56:00 






2004-02-26 




11:16:41 


135 



tion (SED) to the structure of ionized absorbers; and the potential 
for feedback into host galaxies by outflowing nuclear winds. 

The paper is organised as follows. In Section |2] the multi- 
wavelength data are presented, to be followed by a plasma diag- 
nostic approach to the line analysis in Section |3] In Section |4] we 
combine our multi-wavelength data with additional ISO and IRAS 
archived observations to produce an observed SED that is used to 
determine a theoretically motivated ionising spectrum affecting the 
warm absorber properties; this then is used to generate ion popu- 
lations with XSTAR for our spectral fitting. This will set the stage 
for Section |5] considerations on the warm absorber behaviour as 
established by thermodynamic stability arguments (Section 15. It , 
and complex line-of-sight geometry through dust (Section [5.2t . We 
adopt Ho = 73 km s"^ Mpc"\ VIm = 0.25, and = 0.75 
throughout this paper. 



2 OBSERVATIONS & DATA REDUCTION 

As part of a multi-wavelength campaign to better understand the 
global physical processes, absorber kinematics, and geometry of 
IRAS 13349, we observed this nearby quasar using spectrographs 
on Chandra (X-ray), the Hubble Space Telescope (ultraviolet), the 
Hobby-Eberly Teles cope 8-meter (optic al), and the Spitzer Space 
Telescope (infrared; I Werner et al .1120041) . Observations performed 
with Chandra , HST and HET were nearly simultaneous and of 
comparable spectral resolution {R ~ 1000), while "high" (R ^ 
600) and low-resolution (R ~ 70 — 120) IR spectra from Spitzer 
were obtained ~ 1.25 years later as part of the GTO program of 
George Rieke. 



2.1 X-ray Observations with the Chandra HETGS 

IRAS 13349 was observed with the Chan dra HETGS (High En- 
ergy Transmission Gr ating Spectrometer: ICanizares et al.l |2005| ; 
IWeisskopf erai]|2002h over several days in 2004 February for a 
total of ~ 295 ks of usable data as summarised in Table [T] Plus and 
minus first order (±1) spectra were extracted using the latest CIAO 
release (CIAO 4.3 with CALDB 4.4.3), starting with the LI (raw un- 
filtered event) files, which we reprocess to remove hot pixels and 
afterglow events. In order to maximise signal-to-noise (S/N), we 
combine plus and minus 1st orders for both the HEG and MEG. 
The resolving power of the HETGS is E/AE = 1000. We focus 
on the ~ 1.5 — 26 A (0.45 — 8 keV) spectral region in this paper, 
with particular emphasis on the soft (0.5-5 keV; 2.5-26A) band of 
the warm absorber. Analysis of the Chandra line spectra was done 
within the ISI^ analysis package. 

1 http://space.mit.edu/CXC/ISIS/ faouck & Denicolall2000t) 
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Table 2. STIS Observations 



UdLd. Ocl IMdlllc 


Grating 


Date 


Start Time 
(GMT) 


Exposure Time 
(s) 


oSwpOlOlO 


G140L 


2004-02-22 


08:09:48 


2000 


o8wp01020 


G140L 


2004-02-22 


09:31:42 


2630 


o8wp01030 


G140L 


2004-02-22 


11:08:22 


2630 


o8wp01040 


G230L 


2004-02-22 


12:44:03 


2630 



2.2 Simultaneous UV Observations witli the Hubble Space 
Telescope STIS MAMA 

Simultaneous HST (Hubble Space Telescope) UV spectra were 
obtained using the Spa ce Telescope Imaging Spectrograph (STIS) 
jWoodgateetalll 19981) . The faintness of IRAS 13349 in the UV 
made the STIS echelle mode an impractical choice, so we observed 
IRAS 13349 through the 52 x 0.5 "slit with the UV MAMA de- 
tectors and the low resolution UV gratings (R^^^ 1000) G140L and 
G230L in a 4-orbit visit on 2004 Feb 22 starting at GMT 08:09:48. 
Table |2] summarises the STIS observations. 

Three exposures with G140L totaling 7260 s and one exposure 
with G230L totaling 2630 s gave a spectrum covering the 1150- 
3180A spectral range. The G140L spectmm has a resolving power 
that varies with wavelength of 960-1440 over its full spectral range; 
for G230L the resolving power is slightly lower, 530-1040. 

We used the extracted one-dimensional spectra as produced 
by the STScI pipeline for our analysis. Our data required only 
a slight zero-point adjustment (< 1 pixel, corresponding to < 
100 kms~^ in each exposure) to the wavelength scale which we 
determined using the Galactic absorption lines in ou r spectra. Ob- 
servat ions of Galactic H I along this sight-line by iMurphv et al.l 
( 1 19961) show an optically thin column density of 1.07 x 10 cm 
at a mean heliocentric velocity of —19 km s~^. In the region of 
wavelength overlap between the G140L and the G230L spectra, 
the flux levels agree to <0.7%. We have made no adjustments to 
either flux s cale, given that a n absolute accuracy of better than 4% 
is expected l lBostroemll2010l) . 



2.3 Near Simultaneous Optical Observations with the 
Hobby-Eberly 8-meter 

Since IRAS 13349 shows large-amplitude variability at X-ray 
wavelengths, we also obtained a near-simultaneou s spectrum of 
it wit h the 8-meter Hobby-Eberly Telescope (HET; iRamsev et al.l 
Il998l) to constrain its optical properties during the Chandra and 
HST observations. (We note that while the HElljis officially desig- 
nated a 9.2-meter telescope, we conservatively refer to it here as an 
8-meter class given that t his is the average equi valent aperture for a 
typical observation - see lSchneider et al.l20 00l.) HET observations 
were obtained on 2004 Feb 26 w ith the Marcario Low-Resolution 
Spectrograph faill et al.l [r998allbl) . The observations were taken 
with a 600 line mm~^ grating, a slit width of 1", and a GG385 
blocking filter, resulting in a resolving power of 1300 over the ob- 
served wavelength range 4300-7300 A. The total exposure time of 
1283.5 s was split among four sub-exposures. For flux calibration, 
we observed the white-dwarf standard Feige 34 during the middle 
of the night, two hours prior to our observation of IRAS 13349. 



Exposures on a Cd-Ne arc lamp were used to determine the wave- 
length scale, with zero points adjusted using the night-sky lines in 
our spectra. Final flux-calibrated spectra were extracted with the 
standard IRAF reduction codes for single-slit data. 



2.4 Infrared observations with the Spitzer IRS 

Low-resolution (i? = 70 - 120) and high-resolution (R ~ 600) 
infrared spectra were obtained as part of the GTO program of Rieke 
(PIP: 61) on 200 5 June 7 with the Infrared Spectr ograph (IRS, 
Houck etalj|2004h on the Spitzer Space Telescope l lWemer et al.l 
2004h using the standard staring mode. In each spectroscopic mod- 
ule, the source was observed at two nodding positions, and expo- 
sure times were set to 2 x 5 x 14 s in SL (Short-Low) and LL 
(Long-Low), as well as 2 x 5 x 120 s and 2 x 10 x 60 s in SH 
(Short-High) and LH (Low-High), respectively. 

The low-resolution spectral data (LRS) were reduced, from 
the basic calibration data (BCD) files, following the standard pro- 
cedure given in the IRS Data Handboofl The basic steps were 
bad-pixel correction with IRSCLEAN v2. 1, background subtraction 
using the off-source slit, as well as optimal extraction and wave- 
length/flux calibration with the Spectroscopi c Modeling Analysi s 
and Reduction Tool suite (SMART v8.2.L lHigdon et alJlioM) , 
for each nodding position. Spectra were then nod-averaged to im- 
prove the signal-to-noise ratio (S/N). 

The flux calibration appeared consistent between the differ- 
ent modules, with differences less than 5%, and this allowed us 
to combine all the spectra to cover the 5.2—35 fim spectral range 
(4.7—32 iim rest frame). Nevertheless, the continuum at 25 nm 
showed a multiplicative offset from the IRAS 13349 photometric 
flux at 25 ^m obtained from the IRAS Point Source Catalog 
(hereafter, IRAS-PSC). Unfortunately, there is no associated MIPS 
24 fim photometry during the epoch of the Spitzer IRS observa- 
tion, so we cannot rule out the possibility that the object has varied 
in this energy band. However, given that IRAS 13349 is radio quiet, 
and that the shape of the SED suggests that the mid-IR is produced 
by themal emission from dust heated by the central engine, we 
assume that the mid-IR has not varied. We therefore scaled the en- 
tire spectrum by a factor of 1.077 to match the (un-color-corrected) 
IRAS-FSC flux density measurement of 0.84±0.06 Jy at 25 /im. 

We reduce the high-resolution spectral data (HRS) in much the 
same way as previously described for the LRS, with the exception 
that full-aperture instead of optimal extraction is employed, the dis- 
tinction being that pixels along the dispersion axis have equal (full- 
aperture) versus varied (optimal) weighting. Moreover, we could 
not remove background by on-source/off-source subtraction, due 
to the lack of off-source HRS background files in the IRS archives 
associated with the epoch of this observation. Therefore, we extract 
SLl, LL2, and LLl backgrounds from the low-resolution data and 
fit with first-order polynomials. We then use these fits to remove the 
background contribution from the high-resolution spectra, giving a 
good match between the SH and LH continuum flux density level. 
Finally, the overall high-resolution spectrum is scaled by 1.077 to 
match the continuum flux density level at 25 nm to the IRAS 13349 
IRAS-PSC photometric flux density at the same wavelength. 



http ://w w w. as .utexas . edu/mcdonald/het/he t_ gen.0 1 .html 



^ http://ssc.spitzer.caltech.edu/irs/irsinstruiTienthandbook/IRS_Instrament_Handbook.pdf 
* http://tdc-www.harvard.edu/catalogs/iras.html 
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Figure 1. Curve-of-growth calculations representing n = 2 and n = 9 transitions for H-like Is-lsnp and He-like Is^-lsnp transitions for different values 
of the turbulent velocity width h. Plotted data point represent A'j values measured based on a simultaneous fit to all the lines detected for a given H- or He-Hke 
resonance series, assuming h = 100 km s^^ . E.g. plotted points represent a fit (with oscillator strengths locked in the appropriate ratios) to the detected 
n = 2 to n = 9 lines of He-like Is^-lsnp O VII whereas a fit to the H-like Mg XII only includes a fit to the n = 2 and n = 3 lines. (See Section lTll and 
Table[3]for details.) The purpose of the exercise is to demonstrate the importance of the higher order lines for establishing true A'j derivations, where fits to 
single fines (especially the n = 2 a-transition which may lie on the saturated part of the COG) may give an A'j with possible values that span up to two orders 
of magnitude. 
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Table 3. Ionic columns densities based on fits to Chandra observed lines from the full transition series 



Ion 


'^detected 


Vbluc 

(km s"-"^) 


b 

(km s^^) 


logNj 
(cm~2) 


= b 
(km s"-*-) 


i logNj 
(cm-2) 


3 Nh 

(10^1 cm-2) 


O VII 


2 - OO 


_7-„ + 103 
'^''-126 


1371?^ 


18.1 ± 0.5 


100 


1 7 ( 
1 / .f 


^ i 0.5 


1.48 


VIII 


2 - oo 


O»3_107 


288l«» 




100 


17. 


c+0.4 
■6-0.5 


0.94 


Ne IX 


2 - OO 


700 + 131 
-'»»_168 




17 0+0.1 


100 


17. 


7+0.7 
'-0.5 


8.15 


Ne X 


2-6 


000+354 


284ti«^ 


17 + 0.3 

^'■■'-0.4 


100 


17. 


4+0.4 


4.08 


Mg XI 


2-4 


-578^212 


10U_j4g 


1 7 4+3.8 


100 


17. 


4+0.4 
'^-0.5 


13.0 


Mg XII 


2-4 




133t?^ 


17 4+0-4 


100 


17. 


4+0.4 


13.0 


Si XIII 


2-3 


cro7+156 


48l^f 


+ 


100 


17. 


g+1.0 
"-1.1 


22 



See corresponding curve-of-growth plots in Figure [T] 

" Transitions in series which contribute to the best fit TVj ; 00 coiTespond t o series limit . 

^ Outflowing velocity measured against the laboratory wavelengths tabulated bv lVemer et al.l JT996I) 
Turbulent velocity width (km s~^) as determined from the series fit. See il3.1l for discussion on caveats. 
Corresponding ionic column density (cm"^). 
Turbulent velocity width frozen at 100 km ^ . See i|3.1l for discussion on reasoning behind this. 

^ CoiTesponding ionic column density (cm~'^). 

3 Corresponding equivalent Hydrogen column assuming the solar abundance values given bv lWilms etai] l |2000l ). Note however that the numbers assume an 
ionization fraction for all ions to be 50%; as such, quoted numbers are only rough guides, and not to be directly taken at face value. 



3 DATA ANALYSIS AND LINE DETECTIONS 

The intent of this paper is to provide the highest spectral resolution 
multi-wavelength characterisation of the IRAS 13349 spectrum to 
date, as a step toward a more detailed understanding of the structure 
of this quasar and its host galaxy. For this paper, we employ both a 
plasma-diagnostic approach to the anal ysis (this section), as well as 
photoionization modeling with XSTAR jKallman & Bautistall200lh 
for a global analysis (<j4l( of the spectra based on the observed multi- 
wavelength spectral energy distribution (SED). 



3.1 The X-ray Chandra HETGS spectrum 

We begin initially with a plasma diagnostic approach to the line 
fits. To search for individual absorption and emission lines, HETGS 
spectra were binned by two, corresponding to half the spectral res- 
olution, respectively 0.006A and 0.012A FWHM, for the HEG and 
MEG. For these fits, we use a fifth order polynomial to fit local 
continua, rather than a phenomenological broad-band continuum. 
Such an approach allows us to remove the broader fluctuations, 
thereby maximising our ability to measure the narrow lines with- 
out assumptions about details of the continuum, which can be com- 
plex. Using this analysis procedure, we find the X-ray spectrum to 
be dominated by an ionised absorber giving rise to prominent H- 
and He-like absorption lines of oxygen, neon, magnesium, and sil- 
icon, as well as absoiption lines from a variety of iron ions, at a 
bulk outflow velocity in the range —900 kms^^ < I'x-ray ^ 
-500 kms"^ . 

Since a comprehensive photoionization treatment of the 
IRAS 13349 absorber as tied to its SED and theory considerations, 
is explored in detail in Section |4l we concentrate here primarily 
on the key H- and He-like ions as a model-independent way of ac- 
cessing the QSO hot plasma conditions for comparison with later 



analysis. For many of these key ions, high order transitions, corre- 
sponding to a high column density warm absorber is detected (see 
Table[3}. For these key ions, we derive the individual ionic column 
densities (A'j ) based on a fit to the entire detected resonance series 
of lines from n=2-xx, where xx is the highest transition detected; 
e.g. H-Iike Is — np Mg XII and He-like Is^ — IsnpMg XI are 
detected to n = 4 so that our fits to those series are to the Ly/He- 
(ct, P, 7, S) lines with thei r oscillator strength s locked to their tabu- 
lated value as found e.g. in lVemer et al. The individual lines 
in the series are fit using Voigt profiles, although the derived A'j is 
based on a simultaneous fit to all the detected lines in the series. In 
this way, we ensure the best possible determination of A'j (Table[3]l 
by reducing the possibility for false estimates resulting from satu- 
rated or contaminated lines. (Both saturated and contaminated lines 
can falsely indicate a line width and strength which is broader and 
has higher flux than what it truly is for the ion, thereby resulting 
in lower estimates for any given ionic column.) We note however, 
that despite these precautions, the spectral resolution of the Chan- 
dra HETGS, while the best available to date, is still not sufficient 
for probing at the level of measuring the 10 — 20 kms~^ thermal 
widths of these ions. As such, A'j estimates may only be lower lim- 
its. For thoroughness, we adopt a value for the turbulent velocity 
width h — 100 kms~^ to approximate the Chandra HETGS spec- 
tral resolving capabilities, for deriving the ionic columns noted in 
Table [3j we also detail the associated ionic columns based on fits 
where b is allowed as a free parameter. For Ne IX (and to a lesser 
extent Ne X), Fe contaminates many of the stronger lines in the 
series, resulting in a larger measured b, and hence smaller A'j . 

It is also clear that the measured value of A'j strongly depends 
on b as best illustrated by the curve-of-growth (COG) plots in Fig- 
ure [T] Here, we plot COG calc ulations (whic h include the damping 
parameter as presented e.g. in ISpitzeill 19781) for different values of 
the turbulent width, which range from b = 10 - 20 kms"^ (^ 
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Figure 2. Velocity spectra obtained by co-adding ttie 5 strongest resonance lines for H and He like ions of nitrogen, oxygen, neon, magnesium, aluminium, 
silicon, and sulfur. For N, O, Ne, and Mg ions, MEG grating spectra are co-added, whereas HEG data are used for Al, Si, and S. We bin at 200 kms^^ in 
velocity; S is binned at 400 km s^^. - Velocity bins are chosen to reflect the approximate FWHM resolution of the HETGS at the energies where these ions 
are found. To guide the eye, vertical lines denote velocities corresponding to 0, -1000 and -2000 kms^^. Multiple absorption components appear strongest 
for O, Ne, Mg, and Si. The eiTors are Poisson and each bin is independent despite appearances to the contrary. 



the thermal width value of the ion), to 6 = 100 kms~^ (maximal 
resolving power of the H ETGS), to b = 1000 kms^^ (the ap- 
proximate value derived bv lSakoetal.l200ll in their analysis of the 
XMM RGS spectrum). For each ion, two sets of curves correspond- 
ing to the 1 — >■ 2 (Ly / He a) and 1 — !> 9 (representing the series 
limit) transitions are shown. Over plotted as individual points on 
the curves are the values for the different transitions which are de- 
tected in our data, and which contribute in the series fitting of the 
A'j values listed in Table[3] for fixed values of & = 100 kms^^ . 
It is clear that for any given ion, a detection of the higher transition 
series better aids in the determination of A'j . This is well illustrated 
in a comparison of the H- and He-like oxygen series which are de- 
tected to the series limit versus the H- and He-like silicon series, 
where only the a (1 — 2) and /3 (1 — > 3) transitions contribute 
significantly to the fitting. H- and He-like Ne have been excluded 
from the figure due to significant contamination from Fe. For the 
other ions in the figures, contamination to the line series comes 
from other velocity components of the same ions - e.g. for He- like 
O VII , components 1, 2, 3 of the warm absorbers discussed in i ]4.3l 
all contribute to the line series at different velocities. 

To assess further the presence of multiple velocity compo- 
nents, we generate velocity spectra based on the 5 strongest reso- 



nance transitions for the most prominent X-ray ions, namely H- and 
He-like ions of nitrogen, oxygen, neon, magnesium, aluminium, 
silicon and sulphur, (note that not all these are necessarily de- 
tected.) For each individual absorption line, a velocity grid from 
-4000 to 4000 kms^^ is generated centered around the rest wave- 
length of the particular line, i.e. zero velocity corresponds to the 
rest wavelength of the line of choice. The HETGS spectra of counts 
vs wavelength (initially binned by 4) is remapped (through inter- 
polation) to convert to the velocity space. Standard eiTor propaga- 
tion rules are invoked and the remapping retains the Gaussian na- 
ture of the eiTors (i.e. ^/(counts)). The velocity spectra of the five 
strongest resonance lines are then combined (errors are added in 
quadrature) to represent the velocity spectrum of the concerned ion. 
Figure 12] shows the velocity profiles derived in the aforementioned 
way. The velocity is binned at 200 kms^^ except for the bottom 
most panel (i.e. for sulphur) where the binning is at 400 km ^ . At 
wavelengths corresponding to the sulphur lines, which have been 
considered, the velocity resolution for HEG is > 400km s^^. To 
achieve the best results we have used MEG spectra for N, O, Ne 
and Mg while resorting to HEG for Al, Si and S having higher en- 
ergy transitions. The figure clearly shows detection of absorption 
for O, Ne, Mg, Al and Si. The N ions show a hint of absorption, 



© 0000 RAS, MNRAS 000, 000-000 



IRAS 13349+2438 7 



Observed energy (keV) 
6.5 6 5.5 




2 2.2 
Observed wavelength (A) 



2.6 



Figure 3. The Chandra HETGS spectra at > 6 keV reveal two prominent 
emission lines, one due to Fe I Ka at rest energy 6.4 keV and another 
similarly strong line at ~ 7 keV that may be identified with H-like Fe XXVI 
or Fe I K/3 blended with a broadened component of the Fe I Ka line. 



Lya 





^(A) 

Figure 4. HST/STIS low-resolution spectrum of IRAS 13349+2438 in ob- 
served wavelengths. Prominent emission lines are marked. 



whereas no absorption is detected for S ions. We have also looked 
for absorption in carbon, argon and calcium with no significant de- 
tection. For most of the ions the velocity profiles do not have a 
symmetric Gaussian distribution indicating the presence of more 
than one velocity component. We will return to a discussion of the 
X-ray lines in H4.3\ 



3.1.1 Ionized and neutral emission lines 

By far the strongest detected narrow X-ray emission line is the 
Fe I Ka fluorescent line at ~ 6.4 keV (rest) with flux Fkq ~ 
4 X lO"*^ ph cm s (Fig.[3l also Fig. fTSl . Another strong line 
of approximately equal strength appears at rest energy ^ 7 keV 
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Figure 5. Detail of the STIS G140L spectrum of IRAS 13349+2438, in ob- 
served wavelengths, in the Lyo region (top panel) and the C IV region (bot- 
tom panel). The black curves are the obsei'ved data; the red curves are the 
best fit to the spectmm. Emission features ai'e labeled above the spectra, 
and absorption components are marked below. The "G" labels denote fore- 
ground Galactic interstellar absorption lines. 



(Table [§). This energy is best matched to Fe I K/3 but this gives 
an unphysical Ka:K/3 ^ 1, when the ratio shoul d be ~ 10. How- 
ever, given the claim by iLonginotti et al.l (|20o3) based on XMM 
data for a complex iron line that has a narrow component on top 
of a relativistically broadened line, the most likely explanation for 
the measured similarities in the Ka and K/3 line fluxes here is that 
what we have measured in the Chandra data is only a small portion 
of the broadened line. Since a discussion of relativistic effects is 
not the goal of this paper, we will leave it at a reporting of the line 
measurements at high spectral resolution. An alternative possibil- 
ity, especially given the line strength is that the emission is in fact 
due to Fe XXVI Ka at 6.96 keV. We have also checked the ~ 6- 
8 keV spectral region for emission from other ions (e.g. 6.7 keV 
Fe XXV Ka, 7.5 keV Ni I Ka, and 8.26 keV Ni I Ka), but find no 
significant detections. 
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Figure 6. The HET spectram of IRAS 13349, in observed wavelengths, 
shows a typical quasar spectram dominated by broad Balmer emission lines 
and [O in], with the exception of particularly strong Fe II. The lone absorp- 
tion can be attributed to atmospheric Franhofer B-band which we have not 
corrected for, since it has no effect on the conclusions derived from our 
spectram. 

3.2 The UV HST STIS-MAMA Spectra 

To measure the UV fluxes, widths, and redshifts of the emission an d 
absorption lines, we used the IRAlQ task SPECFIT l lKrisslll994| y 
As discussed in i)5.2l the full continuum is a complex mix of di- 
rect light from the active nucleus that is heavily reddened com- 
bined with light scattered from dust, and perhaps free electrons 
also, that also suffers some extinction. We were unable to arrive 
at a model that fully matched the observed shape of the full con- 
tinuum, so, to characterise the emission and absorption lines, we 
used an empirical approach that simply fit a power law that was lo- 
cally optimised around each emission-line blend. Nearly all emis- 
sion lines required both a broad 2500 kms^^ in width) and a 
very broad (> 7500 km s^^ ) component to obtain an adequate fit. 
Table|4]gives the fluxes, velocities, and full-width at half-maximum 
(FWHM) for the fitted emission lines. Velocities are relative to the 
systemic redshift of 2; = 0.10853, based on the observed redshift 
of the 03 A5007 emission line in our HET spectrum (see t]3.3l for 
details). The tabulated line widths are corrected for instrumental 
resolution by subtracting the resolution in quadrature from the mea- 
sured widths. Wavelengths < 1720 A have the resolution of the 
STIS G140L grating as given in the STIS Instrument Handbook, 
ranging from 0.93 A at 1350 A to 0.81 A at 1700 A. Longer wave- 
lengths have the resolution of the G230L grating, which is 3.40 A. 

Figure|4]shows the merged UV spectrum of IRAS 13349 from 
1150-3180 A with the most prominent emission lines labeled. At 
the scale of this figure, Galactic and intrinsic absorption lines are 
not easily visible, aside from Galactic Ly a absorption at 1216 
A, which is blended with geocoronal Ly ol emission. Otherwise, 
our HST spectrum of IRAS 13349 shows two prominent blue- 
shifted absorption-line systems in Ly a, N V, and 4- To measure 

^ IRAF jhttp://rraf.noao.edu7) is distributed by the National Optical As- 
tronomy Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the Na- 
tional Science Foundation. 



Table 4. UV & Optical Emission Lines in IRAS13349-I-2438 



Line 




Flux*" 


Velocity*^ 


FWHM 




(A) 




(kms-i) 


(km s^-^) 




Hubble Space Telescope UV Spectram (Figure|5) 




Ly a. 


1215.67 


5.3 ± 0.1 


-710 ± 368 


3026 ± 150 


Ly a 


1215.67 


10.2 ± 0.1 


-710 ± 368 


1 nni -1— 1 /I A 
10018 ± 140 


N V 


1240.15 


1.1 ± 0.1 


-1130 ± 138 


3189 ± 210 


N V 


1240.15 


2.0 ± 0.3 


-1130 ± 138 


lUUlo ± 14U 


ol IV+U iVJ 


1400 


1.4 ±0.2 


-801 ± 203 


2554 ± obU 


01 IV+U IVJ 


1400 


0.8 ±0.4 


-801 ± 203 


7568 ± 730 


C IV 


1549.05 


2.6 ±0.7 


-1358 ± 105 


2555 ± 360 


C IV 


1549.05 


5.7 ±0.6 


-1358 ± 105 


7568 ± 730 


He II 


1640.70 


0.5 ±0.3 


-1358 ± 105 


2498 ± 360 


He II 


1640.70 


1.8 ±0.6 


— 1358 ± 105 


7549 ± 730 


III] 


1663.48 


0.3 ±0.2 


— 1358 ± 105 


2500 ± 360 


III] 


1663.48 


0.1 ±0.6 


— 1358 ± 105 


7550 ± 730 


Al III 


1857.40 


0.4 ±0.2 


"1 1 K r\r\ 

—891 ± 500 


2174 ± 2190 


Si III] 


1892.03 


0.8 ±0.4 


—210 ± 190 


1482 ± 360 


C III] 


1908.73 


2.4± 1.3 


n 1 r\r\ 

—380 ± 100 


2474 ± 590 


Si IIl]+Cni] 




3.2 ±0.6 




8987 ± 7500 


Mg II 


2798.74 


6.5 ±0.4 


— 318 ± 100 


2130 ± 120 




Hobby-Eberly Spectram 


(Figure |6) 




H e 


3970.07 


2.0 ±0.8 


20 ± 30 


io (D ± 4U 


H € 


3970.07 


2.7 ± 1.2 


-140 ± 30 


5Uo5 ± 4U 


H 


4101.73 


0.5 ±0.5 


20 ± 30 


io/D ± 4(J 


l_i r 


4101.73 


5.4 ±0.5 


-140 ± 30 


5Uo5 ± 4U 


H 7 


4340.46 


0.1 ±0.5 


20 ± 30 


lo7o ± 4(J 


H 7 


4340.46 


8.0 ±0.7 


-140 ±30 


5Uo5 ± 4U 


Fe II 


4434^684 


34.1 ±0.7 


264 ± 24 


2514 ± bo 


He II 


4686.74 


< 0.8 


-140 (fixed) 


io ( u (nxed) 


H fi 


4861.32 


13.4 ±0.6 


20 ±30 


1878 ± 58 


H 13 


4861.32 


21.6 ±0.5 


-140 ±30 


08U 1 ± 142 


[0 III] 


4958.9 


0.41 ± 0.05 


± 30 


421 ± 45 


[0 III] 


4958.9 


2.2 ±0.07 


-330 ± 30 


2556 ± 154 


[0 III] 


5006.8 


1.23 ±0.15 


0±30 


421 ± 30 


[OIII] 


5006.8 


6.4 ±0.21 


-330 ± 30 


2556 ± 154 


[Nil] 


6548.1 


0.6 ±0.1 


0±30 


421 ± 30 


[Nil] 


6548.1 


8.0 ± 1.7 


-330 ± 30 


2650 ± 70 


[Nil] 


6583.4 


1.8 ±0.3 


0±30 


421 ± 30 


[Nil] 


6583.4 


24.0 ±5.1 


-330 ± 30 


2650 ± 70 



We note that the fits include the absorption lines, so there is no fitting bias 
that would artificially blueshift the centroids. Furthermore, the absorption 
fines are too weak to significantly impact fitting results even if unaccounted 
for. "Rest wavelengths ai'e in vacuum for A < 3200 A and in air for A > 
3200 A. ''Observed flux in units of 10"^"* erg cm ^ s ^. '^Velocity is 
relative to a systemic redshift of z = 0.10853, the redshift of the [O III] 
A5007 emission line. 

the equivalent width, position, and FWHM of each line, we use 
simple Gaussian absorption lines in our model to fit the spectrum. 
Figure |5] shows full-resolution plots of the Ly a, N V, Si IV, and 
4 regions of the spectrum overlaid with the best-fit model. Table |5] 
summarises our measurements for each of the detected lines. Line 
widths have been corrected for the instrumental resolution by sub- 
tracting the resolution of 23 1 kms^^ in quadrature from the fitted 
value. We also quote 2a upper limits for the Si IV transitions at 
the same velocities as the other detected components since these 
are useful in constraining the ionisation state of the absorbing gas. 
The lines are slightly broader than the resolution of the L-mode 
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gratings, with intrinsic widths that are consistent with the Doppler 
parameter of 100 kms~^ found for the X-ray absorbers. Assum- 
ing b = 100 kms~^ , we obtain the ionic column densities given 
in the last column of Table[5] Since the intrinsic widths of the UV 
absorption features are not broader than b = 100 kms~^ , these 
column densities can be considered lower limits for Ly a, N V and 
4. The highest velocity component, at — 950kms~^ , has roughly 
the same velocity as the bulk of the X-ray absorption. The lower- 
velocity component, at — 75kms~^ , has a lower outflow velocity 
than most detected X-ray features. 



3.3 The optical Hobby-Eberly 8-in spectrum 

The flux calibrated HET spectrum of IRAS 13349 shown in Fig- 
ure |6]reveals a fairly typical quasar, where broad Balmer and Fe II 
emission lines and narrow 03 emission lines are superposed on a 
blue continuum. Given the luminosity of IRAS 13349, any star- 
burst component or contribution from the host galaxy is completely 
overwhelmed by the QSO emission. 

As we did for the HST spectrum, we use SPECFIT to mea- 
sure the emission lines and continuum in the HET spectrum. For 
the continuum we use a power law in Fx- For the emission lines 
(excluding Fe II ) two Gaussian components are needed to charac- 
terise the line profile. For H/3 and the higher-order Balmer lines, 
we tie the velocities and widths of all the components together 
This is necessary and helpful in deblending these features from the 
ubiquitous Fe II emission. To fit t he complex Fe II emissio n it- 
self, we use the template derived bv lVeron-Cettv et al.l ( |2004|) and 
convolve it with a Gaussian to match with the broadened width 
observed in our spectrum. The best-fit to the continuum yields 
Fa = 3.83 X 10"" x (A/5653. 5)° '^^ ergcm"^ s'^A"^ . Table 
|4]lists the best-fit parameters for the emission lines in the spectrum. 

In addition to these best-fit parameters, we have made some 
empirical measures of the observed H/3 profile and the Fe II 
emission for use in evaluating the mass of the central black hole 
in IRAS 13349 ( i]4.2.3t and for interpreting Eigenvector- 1 (com- 
prised of the relative equivalent widths of 03 A5007 and Fe II 
A4434 — 4684 relative to H/?) as described by iBoroson & GreenI 
( Il992h . To measure the dispersion of the H/? profile directly from 
the data, we subtracted the fitted continuum, the fitted Fe II emis- 
sion, and all other fitted emission lines from the original spec- 
trum. We then computed the empirical full-width at half maximum 
(FWHM) and the dispersion of this net spectrum over the 5266- 
5500 A (4750^962 A, rest) wavelength range of the profile, 
to obtain FWHM= 2468±4.7 kms"^ , velocity dispersion cth/3 = 
1948 ± 3.7 kms"^ , and equivalent width Whib = 88.1 ± 4.5 A. 
(Note that the aforementioned H/3 values are derived using a more 
complex decomposition of the two components, i.e. a narrow core 
which dominates the empirical FWHM, and a broader base needed 
to describe the H/3 line profile (see Table |4] and Figure |6). Sim- 
ilarly, we subtracted the fitted continuum and all other emission 
lines from the original spectrum to obtain the observed Fe II spec- 
trum. Integ rated over the 4434^684 A rest-wavelength range as 
defined by iBoroson & GreenI l ll992l) . we obtain an Fe II flux of 
3.41 X 10"^^ erg cm"^ s"\ and Wpg „ = 82.4 A. For 03 A5007 
we obtain W^q jjjj — 14. OA. In the context of Eigenvector- 1 and 
Eigenvector-2 charact eristics as initia lly discussed by Boroson & 
Green (1992; see also lBorosonll2002h . IRAS 13349 behaves as a 
quasar at one extreme of the Eigenvector- 1 correlations, e.g., the 
ratios O III:H/3 = 0.22, He II:H^ < 0.022, and Fe II:H/3 = 0.97. 
Using the 03 lines in IRAS 13349, we also update the redshift to 
z = 0.10853 ± 0.0001. We note that the very high S/N HET spec- 



trum enables the determination of the centroids of discrete spectral 
features to much better than a resolution element. For the sharp, 
bright [O III] 5007 line, we can determine the centroid to an accu- 
racy of 1.5 pixels corresponding to 30 kms^^ . 

The strong Fe II and weak 03 emission that characterise ex- 
treme Eigenvector- 1 objects are common to Narrow-Line Sey fert 1 
(NLSl) galaxies (e.g. , see JPoggd I2OIII and iBorosonI I2OIII) , and 
IVeron-Cettv & VeronI ( |2006|) classifv IRAS 13349 as a NLSl 
based on its line width. However, our HET spectrum shows that 
IRAS 13349 does not satisfy two out of the three criteria primarily 
used to classify objects as NLSls. The width of H/3 is too broad and 
exceeds the 2000 kms~^ boundary, and the width of H/3 is also 
much broader than that of 03 by nearly a factor of three; NLSls 
should have permitted and broad lines of similar width. Therefore, 
despite its extreme Eigenvector- 1 characteristics, we do not con- 
sider I RAS 13349 to be a NLSl , in agreement with the determina- 
tion of lGrupe et al] ( l2004l20lol) . 

The HET and HST spectra are very similar to the prior obser- 
vations of IRAS 1 3349. - The continuum flux density is nearly iden- 
tical to that observed by IWifls et al.l (Il992h . and ~ 15% brighter 
than that observed bv lHines et alj ( 1200 ih . Our emission-line fluxes 
are also comparable (but note that the units for the last column of 
Table 3 in Hines et al. 2001 should be 10^^* erg cm^'^ s^^), and 
agree to within 15% for H a. However, there are differences of 50- 
100% for lines such as H /J and [O III] , which are badly blended 
with the optical Fe II multiplets. Although we have taken care to 
model explicitly the full Fe II emission spectrum, as have Hines et 
al. (2001), the fits can be highly model dependent, especially when 
one includes multiple components as we have for some lines (such 
as the Balmer lines, [O III], and [N II]). We suspect that the source 
of our differences in line fluxes are due to the methods used in de- 
blending these features. 



3.4 The infrared Spitzer IRS spectra 

As shown in Figure |7] the mid-IR continuum rises fairly smoothly 
to a peak of 0.8 Jy at 30 fim. Superimposed are weak, broad spectral 
features from silicate emission at 10 fim and 18 fim. Weak poly- 
cyclic aromatic hydrocarbon (PAH) resonances are also detected at 
7.7 fim and 1 1.3 f im, indicating a weak starburst contribu tion to the 
mid-IR emission jSchweitzer et alj2006l ; IShi et al.ll2007l) . 

3.4.1 The Low Resolution Spectrum 

We fit the low -resolution spectm m using the IDL_^-fitting rou- 
tines PAHFIT l lSmith et aljr2007i) and MPFIT jMarkwardt.2009,fl 
modified to include a silicate emission component. The continuum 
model consists of 6 fixed-temperature blackbody dust emission 
components (T = 25, 50, 100, 200, 300, 400 K) and a blackbody 
stellar component (5000 K). The silicate emis sion is modeled as 
the product of a Galactic silicate opacity curve jSmith et al.ll2007t) 
and a single-temperature blackbody. This assumes that the silicate 
emission features come primarily from an optically thin region. The 
flux ratio of the 10 /im-to-18 fim silicate features depends strongly 
on temperature, as does the shape of the 10 /im feature. The best-fit 
temperature for the silicate emission is ~ 200 K. Adding additional 
temperature components did not improve the fit. 

The silicate emission temperature is considerably lower than 
the sublimation temperature (^^ 1400 K). This is not surprising 
since the silicate emission features are emitted preferentially over 
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Table 5. UV Absorption Lines in IRAS 13349+2438 



Feature 


Comp 


■^vac 




Av" 


FWHM 






# 


(A) 


(A) 


(km s^i) 


(km s"^) 


(cm '^) 


Ly a. 


1 


1215.67 


0.59 lb 0.05 


- 72 ± 112 


204 ± 19 


14.20+?V?« 

— U.Uu 


Ly ol 


2 


1215.67 


0.31 lb 0.11 


-965 ± 25 


270 ± 19 


13 83+^ °? 


N V 


J 








970 ^ 1 Q 


i4.uy_Q 22 


IN V 


1 


1 Oyi T on 
iZ4Z.?SU 


nil 1 n 1 n 
U.li it U.IU 


— ioy it oUo 


z / u ± ly 


14.04_;^ 


N V 


2 




U.49 ± U.15 


—935 ± 111 


270 ± 19 


14-49„o20 


N V 


2 


1 1 /I T on 


U.4i ± U.zU 


—935 ± 111 


27U ± 19 


1/1 -70+0-24 


ol IV 


1 


1 jyj. /o 


< U.ID 


— / z (nxeclj 


11 yj (nxeu) 


< lo.oo 


Si IV 


1 


1402.77 


< 0.16 


-72 (fixed) 


270 (fixed) 


< 13.65 


Si IV 


2 


1393.76 


< 0.12 


-965 (fixed) 


270 (fixed) 


< 13.21 


Si IV 


2 


1402.77 


< 0.12 


-965 (fixed) 


270 (fixed) 


< 13.51 


C IV 


1 


1548.19 


0.44 ± 0.24 


-76 ± 88 


309 ± 19 


i4.i2i[;:l 


C IV 


1 


1550.77 


0.22 ±0.12 


-76 ± 88 


309 ± 19 




C IV 


2 


1548.19 


0.84 ±0.17 


-930 ± 35 


309 ± 19 


14 K1+0.14 


C IV 


2 


1550.77 


0.55 ±0.14 


-930 ± 35 


309 ± 19 


i4.ao_(, jg 



"Velocity is relative to a systemic redshift of 2: = 0.10853. 
''Ionic column density assuming a Doppler parameter of b = 100 km . 



Table 6. Mid-IR emission lines in IRAS 13349+2438 



Feature" 


\ b 

/^rcst 


^vac 




FWHM^ 


Flux-' 










( X 10^^ Jim) 


(xlO^^^ Wcm^^l 


PAH(+H I?) 


7.678±0.007 


7.700 


4.91±0.99 


15.27±0.10 


11.87±2.40 


[S IV] 


10.506±0.001 


10.511 


2.35±0.14 


4.61±0.73 


4.98±0.23 


PAH(+H I+He II+Mg2Si04?) 


11.243±0.001 


11.250 


5.14±0.95 


16.14±0.09 


7.97±1.41 


PAH(+H I?) 


12.686±0.003 


12.700 


0.99±0.31 


2.47±0.26 


1.17±0.37 


[Ne II] 


12.819±0.001 


12.813 


3.15±0.40 


3.10±0.14 


3.51±0.48 


PAH(+H I?) 


14.241±0.006 


14.250 


3.24±1.75 


4.62±0.82 


1.04±0.74 


[Ne v] 


14.323±0.003 


14.322 


2.96±0.53 


7.29±0.39 


3.30±0.52 


[Ne III] 


15.560±0.002 


15.555 


3.39±0.61 


1.96±0.13 


2.48±0.54 


[S III] 


18.726±0.001 


18.713 


5.08±0.46 


4.94±0.14 


3.73±0.32 


[Ne v] 


24.317±0.009 


24.318 


4.65±1.86 


8.30±1.02 


2.23±0.80 


[O IV] 


25.906±0.002 


25.890 


28.35±3.10 


5.92±0.23 


9.43±1.16 




"Feature's name 








''Measured wavelengths (Rest frame) 








'^Laboratory wavelengths (Vacuum) 







''Equivalent widths 
'^Full-width at half-maximum 
^Flux 



a temperature range of ~ 100 — 200 K. The presence of dust at 
temperatures up to (and possibly exceeding) 400 K is indicated by 
the blackbody components of our continuum model. Keck optical 
interferometry at A"-band s uggests near-IR emiss ion from an inner 
radius of 0.92 ± 0.06 pc ( Kishimoto et alj[200Sl) . reasonably con- 
sistent with dust at temperatures near sublimation. The peak of the 
SED (in vFi,) occurs at ~ 20/im, corresponding to a temperature 
of ~ 200 K. However, with a clumpy dust distribution, the bulk of 
this dust emission may still come from a region close to the subli- 
mation radius. 



3.4.2 The High Resolution Spectrum 

Dust signatures include a strong mid-IR excess and emission fea- 
tures pointing to PAHs. However, it should be noted that while we 
have identified the 11.2/im feature with PAH, additional contri- 
butions from H and/or Mg-rich olivines (e.g. forsterite, Mg2Si04; 



see iMarkwick-Kemper et al]|2007t) are also a possibility. Detailed 
modeling to determine the relative contributions of PAH to Mg-rich 
olivines (if they exist in IRAS 13349) is beyond the scope of this 
paper, however. 

In addition to dust features, low-equivalent width, narrow for- 
bidden emission lines from a range of ionisation states of O, Ne, 
and S are detected (see Table |6] and Figure [8). In particular, the 
[NeV] and [OIV] lines originate from gas that is highly photoion- 
ized by the UV continuum from the AGN. Moreover, the [Ne III], 
[Ne 11], and [S III] lines are likely dominated by AGN emission, 
although there may be contribution from star-forming regions. 

iTommasin et al.l ( l201(]h presented a comprehensive mid- 
IR high-resolution spectroscopic survey of 91 Seyfert galax- 
ies and derived useful observational and semi-analytical diag- 
nostics that we also use here to assess the degree to which 
IRAS 13349 is AGN- or starburst-dominated. In particular, a 
good tracer of AGN contribution can be found in a compar- 
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Figure 7. The best fit model (green) to the low resolution Spitzer IRS spec- 
trum (x 1.077) showing all contributions (black body continua, silicate, 
PAH, starlight, and narrow-emission lines). To contrast, grey curve shows 
contribution from black bodies, starlight, and silicate; additional contribu- 
tion from PAH emission (blue), and naiTow-emission lines [S IV], [Ne 11], 
[Ne III] and [O IV] (puiple) are also shown. Model components contribut- 
ing to the aforementioned include silicate emission (bottom, black) super- 
imposed on a black body continuum (red) and hot dust or starlight (magenta; 
note that it is not possible to distinguish the two without observations at 
near-lR wavelengths). The best fit silicate emission temperature is ~ 160 K 
while the continuum components have temperature range 50 — 400 K. 
(The T=25K continuum component has negligible contribution to the fit 
and hence not shown as part of the model here.) 

ison of the [Ne v]14.32/im : [Ne Il]12.81^tm ratio against 
[O iv]24.32^m : [Ne Il]12.81/im. - For IRAS 13349, these val- 
ues are respectively 0.800 ± 0.223 and 2.958 ± 0.322, pointing 
to an 80-90% AGN contribution to the IRAS 13349 IR emission, 
according to Fig. 8a of Tommasin et al. 



4 THE IRAS 13349+2438 SPECTRAL ENERGY 
DISTRIBUTION AND ITS INFLUENCES ON THE 
WARM ABSORBER: OBSERVATIONS AND THEORY 

In having considered the line detections separately in the different 
wavebands, here we assess the IRAS 13349 warm absorber based 
on photoionisation modeling as tied to this QSO's SED. In par- 
ticular, we investigate the impact of different SED components (in- 
cluding the X-ray Compton power-law component, soft excess, and 
big blue bump/accretion disk) spanning UV to X-rays on affect- 
ing warm absorber thermodynamic stability conclusions. We also 
investigate the methods for assessing the IRAS 13349 black hole 
mass and accretion rate based on the shape of the SED. 

4.1 The Observed Spectral Energy Distribution of IRAS 
13349-1-2438 

To assess photoionization scenarios for the absorbing gas in 
IRAS 13349, we first need to know the intrinsic spectral en- 
ergy distribution (SED) of the active nucleus that is illuminat- 
ing the surrounding gas and dust. To start, we assemble an IR to 
X-ray SED based upon our observations and archival data (see 
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Figure 8. High resolution Spitzer IRS LH (bottom panel) and SH (top 
panel) spectra of IRAS 13349. The spectrum shows dust features attributed 
to lO^m (broad continuum bump from 9 — 13/^m) and 18/xm emission from 
dust in silicates. The 18/^m is so broad that it is difficult to distinguish from 
a blackbody peaking at 200 K. 

Figure In t he IR, we start with 170 far IR ISOPHOT 
measurements jSpinoglio et al. I l2002h . while mid-IR (12, 25 and 
60 Aim) IRAS data a re taken from the Faint Source Catalog 
of llMoshir et"aLl ( Il990h . and near-IR photometric po ints are the 
2MASS 1-2 data from the Large Galaxy Atlas of Ijarrett et al.l 
lj2003|). For comparison, our Spitzer IRS spectrum is overlaid on 
these data after scaling to match the 12 /im and 25 fim IRAS pho- 
tometry (as noted in m.4\ . In the optical and UV bands, we correct 
our observed HET and HST spectra for foreground Galactic extinc- 
tion of E{B ~V) — 0.013, which we obtain from NED, assuming 
a mean Galactic extinction curve (fCardelli et al. 1989) with a ratio 
of total to selective extinction of Rv ~ 3.1. Our Chandra spec- 
trum in the X-ray is corrected for foreground Galactic absorption 
as described in §2. 

For comparison, we generate a "generic" composite optical- 
UV QSO spectrum based on t he Sloan Digital Sky Surve y (SDSS) 
composite quasar spectrum of IVanden Berk et al.l ( 1200 iT) (for A > 
3200 A rest) and the ra dio-quiet HST composite quasar spec- 
trum of iTelfer et af] ( l2002h at shorter wavelengths. If we normalise 
this composite spectrum to match our extinction-corrected HET 
spectrum at 5700 A rest (resulting in the solid blue curve of Fig- 
ure|9l hereafter "5700A -normalised-generic-composite") what im- 
mediately stands out is the large deficit of UV and far UV flux 
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in IRAS 13349, which is likely reradiated by heated dust in the 
IR since the observed SED peaks at 30 /im, rest. In the optical 
band, the HET spectrum of IRAS 13349 is virtually identical to 
that of the SDSS composite portion of the 5700A -normalised- 
generic-composite, indicating that the IRAS 13349 optical spec- 
trum is only mo destly reddened or extincted at these wavelengths. 
(The analysis o f lHopkinslliooi shows that the full SDSS compos- 
ite has an internal SMC-like extinction of E{B - V) ~ 0.013 
at most.) The Balmer decrement measured in our HET data is 
2.85 ± 0.18 for the ratio of HQ-to-H/3. Correcting this for fore- 
ground extinction, it becomes 2.817 ± 0.18. To coixect this to a 
nominal Case B value of 2.76 implies E{B - V) = 0.021 for 
an SMC-like extinction law. Comparing the Galactic-extinction- 
coiTected spectrum of IRAS 13349 to the SDSS composite in more 
detail, the ratio (in the continuum) is flat from 4000-6000 A; if 
anything, IRAS 13349 appears a bit bluer. So, internal extinction 
on the order of 0.01 to 0.02 (with an SMC-like law) is consis- 
tent with E{B — V)= 0.021 as derived from the Balmer decre- 
ment. - For this value Ay = E(B - V)_R v = 0.09 ± 0.02 taking 
Rv = 4.05 ± 0.8 ala lCalzetti etal] ( l2000h . 

Based on the 5700A -generic-normalised-composite spectrum 
we find that the UV/optical peak (i.e. the peak of the solid blue 
SED) is 0.6 dex below the peak radiation in IR (at ~ 10'"^ A). 
According t o the median radio-quiet QSO SED in Figure 12 of 
lElvij ( Il994h . the UV peak for a "normal QSO" should be 0.25 dex 
above the IR peak. Assuming IRAS 13349 to be in this category 
(IR and optical spectra suggest this to be the case for IRAS 13349) 
its UV/optical peak is then 0.85 dex below that of an unabsorbed 
"normal QSO". This suggests that either we are viewing the op- 
tical spectrum of IRAS 13349 through a gray screen, or largely 
in scattered light, but from a scattering region that again is largely 
colour neutral at wavelengths longward of 4000 A. At shorter wave- 
lengths, the observed UV spectrum is likely some complex mix of 
scattered and absorbed light from the central regions. The complex 
geometry of this radiative transfer is difficult to unra vel even us- 
ing th e added benefit of polarisation information (W92: lHines et al.l 
l200ll) . so to produce an SED that is fully corrected for these ef- 
fects, we will assume that scaling the composite QSO spectr um 
to the level implied by the median QSO SED in lElvisI i ll 9941) is 
a good representation of the intrinsic SED of the active nucleus 
of IRAS 13349. The resultant SED (see Fig. [9l solid green and 
associated dotted blue line connecting it to the IR and X-rays) 
looks like a qualitatively good description for the intrinsic SED 
of IRAS 13349 since the extreme ultraviolet portion of the com- 
posite spectrum is declining smoothly to higher frequency in a way 
that would provide a good match to the soft X-ray portion of the 
spectrum observed with Chandra . Assuming this to be the in- 
trinsic IRAS 13349 spectrum, at its coiTected UV luminosity of 
£-2500 ~ 8.1 X lO'^^ ergs Hz~^, we derived aox ~ —1.83, a 
value that is —0. 25 lower than the no minal Uox = —1.58 value ex- 
pected from the lYoung et al.l 1 I2OIOI) relation. This slightly lower 
aox we observe for IRAS 13349 is consistent with the weak 4 
absorption that we see in the HST spectrum (c.f. iGallagher et al.l 
l2006h . 

We have also checked the Elvis et al. SED against more re- 
cent compilations of QSO SEDs ( IShand201 iMRichards et alj2006l : 
iHopkins et all |2007|) built from various combinations of SDSS, 
Spitzer, FUSE, and HST data. The variation in the UV-IR peak dif- 
ference (in log(i^F^)) is less than a factor of 2 when comparing 
these different samples. (Note that the Hopkins et al. SEDs use the 
Richards et al. sample of objects but adds an additional Qox correc- 



tion based on the object's peak UV flux; also, the Richards sample 
of QSO, in using primarily SDSS have a median redshift 2 > 0.5). 

To assess luminosity contributions for the different wave band 
components of these SEDs, we integrate the line-segment SEDs 
of Figure |9] over the energy bands noted in Table |7] we assume 
Ho = 73 km s~^ Mpc-\ f^M = 0.25, and f^A = 0.75. 
These luminosities are then used as relevant for calculations in 
subsequent sections, either in whole or in part. The bolometric 
luminosity we derived from the piecemeal integration is Lboi ~ 
3.12 X 10"^ ergs"^ . 



4.2 Reconstruction of the UV to X-ray SED and its 
components: tlieory and observations 

The radiation from the central regions of AGN is likely to peak in 
the extreme ultraviolet at ~ 10 - 100 eV (~ 124 - 1240 A), an 
energy range often suffering from extinction effects due to Galac- 
tic and/or intrinsic dust. As stated in i]4.1| such is the case for 
IRAS 13349, which shows a large UV and far-UV deficit in its 
observed SED, which we attribute to dust extinction, a conclusion 
also supported by Spitzer IR observations (see i]3.4.1l for details). 
Likewise, the two most important spectral components featuring in 
this !(n-observable energy range are the disc blackbody from an ac- 
cretion disc peaking between 10 - 100 eV (typical of AGN), and the 
soft excess, typically manifesting at < 1 keV. Here, using theoreti- 
cal considerations in combination with our multi-waveband obser- 
vations of IRAS 13349, we derive what we believe to be the most 
likely scenario for the intrinsic IRAS 13349 SED, sans extinction. 

The ionising continuum dictating the ionisation state of the 
absorbers in AGN is primarily contained in the UV-to-X-ray range. 
Therefore we focus our discussion on this energy range, with spe- 
cial attention paid to recreating the extincted UV part of the spec- 
trum. 

The general mathematical form we use in building the 
IRAS 13349 theoretical SED is: 

,,3 



2Tvh 



exp{h v/ks Tsi 



+ Adbb fdbbiv, Tij 



+ Adbb fdbb{v, Tin) 



(a) 



(b) 



(1) 

where the first and third terms in both Equations [T]i and[T]) represent 
the 1-10 keV X-ray power-law and the big blue bump (modeled by 
"disk-blackbody") respectively. As can be seen, the two equations 
differ only in the 2nd term, associated with competing models de- 
scribing the soft excess (see i]4.2.1| for details). We proceed with 
a discussion of this component, to be followed by more in-depth 
discussion of the big blue bump: its derivation, and the dependence 
of its shape on the central black hole mass and accretion rate. 

4.2.1 The sub-ke V soft excess 

Observations of most AGN reveal that the SED between 2-10 keV 
is well approximated by a power-law with spectral indices a 0.9 
(photon indices F ^ 1.9). However, if this power-law is extended 
to lower energies (< 1 keV), for most type I AGN, there is addi- 
tional unaccounted radiation, which has come to be known as the 



© 0000 RAS, MNRAS 000, 000-000 



IRAS 13349+2438 13 



10^ 



Kes, (log A) 

10* 1000 



100 



10 



10- 



10- 



o 

lao 
t-, 

(0 

DO 

o 



10- 



10- 



-1 1 1 1 I I 



-1 1 I I I I I — I 1 1 1 I I I I I — I 1 1 [T 



-1 1 1 I I I I — I — I 1 1 1 I I I I — I — I r 



Composite QSO 




^rest (log Hz) 

Figure 9. Observed and intrinsic Spectral Energy Distribution (SED) of IRA S 13349 based on our o bservations. The x axis (both v and A) is in the rest frame. 
Solid points with en'or bars ai'c archived d ata from (1) ISO 17 0/^m data from lSpinoglio et al] J2002| ), as labeled. (2) IR AS 12, 25, 60 ^m measurements of the 
IRAS Faint Source Catalog, Version 2.0 - iMoshir et aljl990l and (3) 2MASS Large Galaxy Adas of ljarrett et alj2003l Spectra are also shown (in solid black) 
of recent observations with the Spitzer IRS (IR), the HET (optical), the HST STIS (UV), and the Chandra HETGS (X-ray). The "Composite QSO" spectrum 
(solid green line) has been renoraialised upw ard by 0.85 d ex so that its EUV peak lies 0.25 dex above the IR peak of the observed IRAS 13349 spectrum, 
as is found in the mean RQ QSO spectrum of lElvisI |l994b. Straight line segments between 10^" Hz and 10^ Hz show a range of simple assumptions for 
what the IRAS 13349 SED can be in the EUV energy range where there is no data. At one extreme, we take our observed SED without any corrections and 
decompose it into a set of power-law segments connecting the observed HST spectrum to the lowest (~0.3 keV, rest) energy Chandra point (solid red). The 
solid blue line represents a "corrected" UV spectrum where the "correction" is derived by assuming that the intrinsic spectrum is represented by the Composite 
QSO spectrum normalised to the observed HET spectrum of IRAS 13349 assuming a break at the peak of the Composite spectrum (1000 A rest) that we again 
extrapolate down to ~0 .3 keV in the X-ray. The dotted blue line follows the continuum established by the Composite spectrum renonnalised by the relative 
UV and IR peaks of the lElvij <1994h mean RQ QSO SED, which we again extrapolate through the EUV using a power law from the peak of the Composite 
spectrum at 1000 A to the Chandra 0.3 keV soft X-ray point. 



soft excess. This component is pa rticularly important for shaping 
the nature of the X-ray absorber dChakravorty et al.1l2012l) . How- 
ever, the physical origin of the soft excess is not understood, and 
presently the community largely treats it as a phenomenological 
spectral component modeled as a blackbody (Equation [T^) with 
temperature Tg^ ~ 100 — 200 eV, which generally does a good 
job of describing most observations. Since this temperature is too 
hot for the temperature, Tin, of the innermost stable circular or- 
bit of the accretion disc in AGN (see Equation |2ll, the soft excess 
is often believed to be a separate spectral component altogether, 
or a reprocessed (shortward in wavelength) extension of the ac- 
cretion disc component. Our fit to the Chandra HETGS spectrum 
of IRAS 13349 in §3.1 gives a temperature for the soft excess of 
Tse = 99 eV. 



soft excess is the thermal Comptonisation model nthcomp 
dLightman & ZdziMskil 1 19871 : IZdziarski et"al1 1 19961: IZvcki et all 
1 19991) included in QxsPEC v.12.5 l lAmaudlll996l)~ In this descrip- 
tion, the seed photons from the accretion disc are reprocessed 
by the thermal plasma to generate sufficient photons at sub-keV 
energies to mimic the soft excess. The high energy cut-off of this 
component is parametrised by the electron temperature Tc, and 
the low energy rollover is dependent on the effective temperature. 
Tin ~ T(Rin) (Equation|2ll of the seed photons from the accretion 
disc. Between the low and high energy rollovers, the shape of 
the spectrum is approximated by an asymptotic power-law with 
photon index Fc ~ fy"^''^, where y is the Compton y-parameter, 
which gives a measure of the extent of Compton reprocessing; 



An alteiTiative model fc{i^) (Equation [TJ)) for the 



http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/ 



© 0000 RAS, MNRAS 000, 000-000 



14 Lee et al. 

Table 7. Luminostiy of IRAS 13349+2438 in different energy bands 



SED 


^ Lbol 


■' LiR 


® Lopt 


h T 

Luv 


' Leuv 


J Lx 




(Xl0""5ergs-i) 


(Xl0*^ergs-i) 


(xlO-'^ergs-i) 


(xlO-l^crgs-i) 


(xlO'^^ergs-i) 


(Xl0'*''crgs-i) 


a 1 


1.42 


1.19 


0.72 


0.36 


0.68 


5.55 




1.74 


1.19 


0.72 


1.44 


2.77 


5.55 


" 3 


3.90 


1.51 


3.72 


9.73 


9.89 


5.55 


d 4 


2.03 




1.95 


9.60 


7.19 


5.42 



"Observed SED - the solid red line in Figure |9] 
''Corrected SED - the solid blue line in Figure|9]and in all panels of Figure [To] 
'^CoiTected SED - the dotted blue line in Figure|9]and in all panels of Figure [Tol 
''Theoretical SED - the solid red line in bottom panel of Figure [Tol 
'^Bolometric luminosity integrated from log v = W.'o — 20. 
^IR luminosity integrated from log v = 11.5 — 14.5. 
f Optical (HET) luminosity integrated from log v = 14.5 — 14.87. 
''Ultraviolet (HST) luminosity integrated from log v = 14.87 — 15.5. 
'Extreme ultraviolet (unobserved) luminosity integrated from log u = 15.5 — 17.0. 
''X-ray luminosity integrated from log u = 17.0 — 20. 



i.e., the larger the value of y, the greater the fraction of photons 
reprocessed from the accretion disc. For our modeling efforts, we 
adopt Tc = Tse ~ 100 eV, consistent with our Chandra soft X-ray 
data for IRAS 13349. The associated normalisation constant Ase 
was determined from the ratio of the 35 — 5 A (0.35 - 2.5 keV) 
photon flux in the soft excess component to that in the power-law, 
as seen in the Chandra HETGS data for IRAS 13349. 

The aforementioned models of the soft excess do not influ- 
ence the UV part of the spectrum and hence do not affect the pre- 
dictions for the shape of the big blue bump (see next section). The 
soft excess is an important component influencing the ions which 
absorb the soft X-ray radiation. Our theoretical simulations show 
that there is a slight difference in the predicted ion fraction of the 
various soft X-ray ions, depending on choice of models for the soft 
excess. However, while the Chandra HETGS data cover the energy 
range of the soft excess, these are not sensitive enough to detect 
these relatively small differences in ion abundances, and hence can- 
not be used to differentiate between the nthcomp or blackbody 
models. Therefore, for the remainder of our analysis we adopt the 
nthcomp model with Tc ~ 2.3 and Tc = 100 eV, acknowledging 
that an alternative blackbody model with Tss = 100 eV is likely 
to produce equally good results. 



4.2.2 The "Big Blue Bump" 

Multi-wavelength observations suggest that AGN continua peak in 
the E UV energy band and usually dominate the quasar luminosity 
(e.g.. lElvisll994h . This spectral component, often referred to as the 
"Big Blue Bump", is considered to be the signature of the presence 
of an accretion disc. Yet, for many systems, it can only be partially 
observed in the UV-EUV energy range. As such, an attempt to re- 
construct it requires careful theoretical considerations . 

According to stand ard thin disc accretion theory 
dShakura & SunvaevI Il973h . radiation from the accretion disc 
may be modeled as a sum of local blackbodies emitted from the 
different annuli of the disc at different radii, the temperature of the 
annulus at radius R being 



T{R) = 6.3 X 10 



5 / 



M 



\mEddJ \W^Mq 
jPetersonll 19971 : [Frank et al.ll2002h where m is the accretion rate of 



the central black hole of mass Mbh, rriEdd is its Eddington ac- 
cretion rate and Rg is the Schwarzschild radius. The normalisation 
constant Adit for this spectral component is given by 



[ Ri„/km 
\d/(10 kpc) 



cos 9 



(3) 



for an observer at a distance D whose line-of-sight makes an an- 
gle i to the normal to the disc plane. Rin is the radial distance of 
the innermost stable annulus of the accretion disc from the black 
hole. Thus, the radiation from the accretion disc has direct depen- 
dence on the mass of the black hole and its accretion rate. As such, 
the shape of the spectral energy distribution can provide important 
diagnostic power for assessing the fundamental parameters of the 
black hole. 

We acknowledge that more rigorous models exist for mod- 
eling the big blue bump, that involve real radiative transfer in 
the accretion disc (se e e.g. iBlaes et al]|200ll : iHubenv et alj|2000l 
200 ll: iHui et alj II2OO5I) and/or the black h ole spin effects (see e.g. 



Davis et al.ll2005l : lDavis & Hubenvll2006h . For the same black hole 



mass and accretion rate, relative to a simple "disk-blackbody", 
these more-involved models change the location of the peak and 
the shape of the EUV spectrum, as UV and EUV flux gets ab- 
sorbed and re-emitted at longer wavelengths in these models. For 
models which include the black hole spin, the peak of the accre- 
tion disk spectrum is pushed to higher energies for increasing black 
hole spin. Thus, qualitatively we can expect that using these models 
would result in higher black hole mass as compared to the best-fit 
results obtained by using "disk-blackbody", for the same accretion 
rate. However, since the UV spectrum of IRAS 13349 is heavily ex- 
tincted, we cannot make observationally based distinctions between 
the different accretion disc models. Given this, in the next section 
we opt for the less-computationally expensive "di sk-blackbody" 
model ( lMitsudaetal]|l984 iMakishima eTall 1 1 9861) to assess the 
mass and accretion rate of the black hole. While we present results 
based on fits to DISKBB to allow for easier comparison to work 
on this and other AGN by other authors, we acknowledge that the 
model EZDISKBB jZimmerman et al.ll2005l) may be the theoreti- 
cally more sound model, and therefore also discuss our subsequent 
SED results based on a comparison of the two, where relevant. In 
brief, the major difference in the two models is most notable at 
R < 107?in, whereby the DISKBB predicted temperature continues 
to increase in contrast to EZDISKBB; this difference stems primarily 
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Figure 10. Theoretical SEDS calculated for IRAS 13349 based on multi-wavelength data shown in Fig. |9] The filled red squares reflect the observed 
IRAS 13349 SED based on ISO, IRAS, Spitz er, HET and Chandra. Gr een squares are the renormalised observed "generic-composi te QSO" spectrum made 
up of a SDSS composite of ~ 2000 quasars l lVanden Berk et aLlliooH) longward of 3200 A and the HST RQ quasar composite of iTelfer et al] l l2002h . The 
soUd and dotted blue lines trace out the same SED as in Fig. |9] (BOTTOM PANEL:) The theoretically computed SED (solid red) that best matches the 
observed IRAS 13349 SED corresponds to a black hole mass of Mbh ~ 10®'^^ Mq, an accretion rate of rii/rriEdd ~ 0-2 and = 10i?s. On the 
other hand, if the accretion disc is completely ignored and one draws a SED based only on the Chandra-HETG spectra, we get the dashed black curve. (TOP 
PANELS:) Theoretical SEDs comparing spectral shape differences for different combinations of Mbh, lii/mEdd and Rin/Ra as discussed in Section I431 
for comparison and edification. 



from the EZDISKBB imposed boundary condition that the viscous 
torque be zero at the inner edge of the disc. 



4.2.3 Determining the mass and accretion rate of the 
IRAS 13349 black hole 

Using the observed IRAS 13349 SE P (see ^43] discussion), in 
combination with Equation [T] and the IShakura & SunvaevI ( Il973h 
model (applied to Equation [T] 3rd term), we generate a series of 
SEDs for different values of ?n and Mbh to obtain the best match 
to the IRAS 13349 SED (Figure[9l gUfor details). We carry out an 
iterative, step-by-step approach to reconstruct the UV-EUV SED, 
with permutations of the relevant parameters, as described below. 

We begin by running models using a coarse grid of 
various combinations of the parameter values rfa/riiEdd = 
(0.05,0.1,0.5,1.0) and Mbh = (10^ 10^ 10^ Mq for each 



rh/riiEdd,) with a fixed Rin = 3Rs. For this initial run, we find 
that rh/riiEdd = 0.05 and Mbh = 9.0 gives a reasonable match to 
the observed IRAS 13349 SED (see top left panel of Figure [Tol 
comparing theoretical models with UV composite QSO spectra, 
green data points). It should be noted that we can expect a range 
of acceptable matches between theory and observations for differ- 
ent permutations of key parameters, primarily due to the fact that 
Equation[T] which we use for modeling the big blue bump, is degen- 
erate between mass and accretion rate. The top panels of Figure [Tol 
demonstrate the effects on the shape of the theoretical SED when 
various combinations of Rin, Mbh, and rii/riiEdd are fixed/varied. 

Therefore, to reduce uncertainties, we use a black hole mass 
derived from the line width (cth/s = 1948 ± 3.7kms~^ ; 
see Section [33t , and the optical continuum luminosity AL5100 = 
6.28 X 10'*'' erg s^*, based on the renormalised, green Compos- 
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Figure 11. Based on grid runs using DISKBB, a tlieoretical representation of al! tlie SEDs that we liave investigated to find the black hole parameters that best 
match the composite QSO spectra normalised to the flux of IRAS 13349. The big blue bump SED can be mathematically described by Tin and A^bt which 
we plot for each combination of the parameters Mbh, rh/niBdd Rin- The three different grids (colour coded and labeled appropriately) of Mbh and 
rii/niEdd correspond to three values of Ri„/Rs. The black star represents the best match SED (red curve, lower panel of Figure llOl for Mbh ~ 10*'^^ Mq, 
an accretion rate of rii/mEdd ~ 0-2 and Ri„ = lORg. See Section l4.2.3l for details. 



ite QSO spectrum in Figure |9] Based on the iMcGill et aP ( l2008h 
relation, 

log Mbh = 7.68+2 log ^ km s'^ +0.518 log ^^^^ erg s-\4) 

we derive for IRAS 13349, log Mbh = 9.10 M©; the McGill 
relation has an rms scatter in log Mbh ~ 0.09. Barring better con- 
straints, we use Mbh ~ lO'"* Mq for input into Equation[T]to de- 
termine the accretion disc component of our theoretical intrinsic 
SED. Given that this Mbh estimate is close to our "best fit" value 
from the initial coarse grid of runs based on a fixed Rin = 3i?s, 
we next run a finer grid of models allowing all the aforementioned 
three parameters free. Figure [TO] (bottom panel) shows the best 
theoretical SED match (red curve) to the UV-EUV flux distribu- 
tions based on our hypothesised extrapolation of observed data in 
the UV-EUV spectral domain (green points). The best match SED 
has optimal parameters Mbh = 10**'^''' Mq, rh/mEdd = 0.2, and 
-Rin ~ 10i?s- Figure [TT] shows Tin versus Adtb for all theoretical 
values of Mbh, rh/riiEdd and Rin that we have investigated, in- 
cluding where in these curves our best fit lies. For thoroughness, 
we also perform the same exercise using EZDISKBB, and find that 
while the "best fit" SED shape remains the same, the same exer- 
cise with EZDISKBB gives a higher accretion rate for fixed val- 
ues of Rin and i\'fBH. From our trials at finding the best match 



based on EZDISKBB, the relevant SED parameters bracket a similar 
best range found using DISKBB, i.e. 10** ™ < Mbh < IO'-'-^'^Mq, 
0.2 < m/riiEdd < 0.7, and 3 < i?in < 10i?s- 



The shape obtained for the big blue bump is then comple- 
mented with X-ray power-law ( i]4.2b . and soft excess (discussed 
in i]4.2. U components to derive the full UV-XRay SED (Figure [Tol 
bottom panel - red) needed for input to photoionization modeling.- 
See luminosity information associated with SED #4 of Table |7] 
(Note that Lir, which is primarily attributed to reradiation by dust 
that has absorbed part of the optical-UV radiation is excluded here, 
since it has no bearings on the actual disc spectrum.) Note that if 
the accretion disc component is completely ignored and a SED is 
constructed on the basis of only the Chandra HETGS X-ray data 
for IRAS 13349, then one would obtain the black dashed curve 
shown in the bottom panel of Figure [Tol Such an SED was used by 
iHolczeret al] ( l2007l) to study the warm absorber in IRAS 13349. In 
Sections r4.2.4| and |5.1| we discuss differences in, and effects on, de- 
rived absorber properties based on an ionising SED that accounts 
for only X-ray contributions versus one which also includes disc 
contributions. 
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Figure 12. Ion fractions versus gas temperature as a function of tlie two 
relevant ionising SEDs. The line styles follow the same scheme used in 
Figure [Tol : the soHd red curves are obtained using the full UV to X-ray SED 
whereas the dashed black curves are a result of an X-ray only SED based on 
only the Chandra-HETG continuum. The panels A-C are for the UV ions 
and D-H are for the X-ray ions as labeled. The dotted vertical lines mnning 
across panels D-H mark the range 2.0 < log^ < 3.0 corresponding to the 
full UV to X-ray ionising SED. This is the range of ^ within which the ion 
fractions for Ne IX, Mg XI and O VIII peak. The importance of this range 
of the ionisation parameter is discussed in Sections [4. Bl and lsTI 



4.2.4 Advantages to using the full UV to X-ray SED for 
assessing the location(s) of the absorbing gas 

Approximately half of all low-z AGN show high-ionisation UV ab- 
sorption lines and signatures of even higher ionisation warm ab- 
sorption in the X-ray data. In assessing the role of the disc ionising 
radiation for affecting our conclusions on UV and X-ray warm ab- 
sorber properties (e.g. common versus distinct origins; continuous 
vs. clumpy clouds), we use XSTAR to generate ion fractions (Fig- 
ure I12t as a function of the temperature of the absorbing gas for 
the UV (panels A-C) and the X-ray (panels D-H) ions, using an 
SED that includes the big blue bump, the sub-keV soft excess and 
the X-ray power-law (Figure [TO]lower panel - red curve) versus an 
SED that only accounts for the X-ray spectral components, namely 
the power-law and the black body soft excess, as determined from 
continuum fits to the X-ray data (Fig.[TO]black dashed curve). 

One of the most important ions in the X-ray warm absorber is 
He-like O VII. As can be seen in Figure[T2l its ion fraction is under- 



predicted by ~ 25% if the big blue bump is excluded from the ion- 
ising continuum (panel D). The relative abundance of O VII and the 
H-like O VIII, another of the important ions in the X-ray warm ab- 
sorber, are sensitively inter- related. Hence under- prediction of the 
O VII ion would be associated with a complementary over pre- 
diction of the O VIII (panel G). Mg XI is also significantly over- 
predicted by a SED which excludes the big blue bump. Further- 
more, it can be seen that the inclusion of the big blue bump pre- 
dicts a larger ion fraction for the UV ions (panels A-C), and hence 
a higher likelihood that they are produced in the same gas as that 
responsible for O VII (panel-D) seen in X-ray absorbers. This is 
contrary to the conclusions one would arrive at, if considering an 
X-ray-only SED. 

4.3 The ionised dusty absorbers in IRAS 13349 

Using XSTAR 2.2.1, we generate "warm absorber" models with 
ion populations that reflect the ionising spectrum of our preferred 
IRAS 13349 SED (see Figure [lOl bottom panel, red SED; ^4231 
for details). We employ the analytic version of XSTAR which enable 
line broadening and optical depth calculations in real time as rele- 
vant to the data. In addition, the granularity associated with analytic 
models are only in ^, which we mitigate by creating population files 
sampled with A(log^) < 0.1 . (In contrast, one can expect table 
models to have granularity in all free parameters. - Private Com- 
munication: Tim Kallman.) 

Using codes we developed for the ISIS jHouck & DenicoQ 

I2OOOI) fitting/analysis package in combination with the aforemen- 
tioned XSTAR-generated models, we initially fit the data (binned 
to the HETGS resolution) with 9 absorber components spanning 
the —4 < log^ < 4 range, and find that it is the high ^ 
(log^ > 1) gas which drive the fit. In also investigating dust- 
less (Section |4. 3. 2| l and dusty (Section [4. 3. 3t absorption, we find 
that a statistically good fit describing the HETGS data (0.5-1.3 
keV MEG and 1.2-8.9 keV HEG chosen to maximise both spectral 
resolution and throughput) is a power-law plus blackbody contin- 
uum absorbed by two ionised absorbers [hereafter WA- 1 (log ^ = 
1.68tj;;lt logNn = 21.031HI) and WA-2 (logC = 3.261^1, 
logNn = 21.66lJ!;JJ^)] at similar (~ 700-800 kms"^) veloc- 
ities, and iron dust (see Section 14.3.31 for details), intrinsic to 
the source. The reduced Cash statistic AC for this model is 
1.147 (2887.127/2540) and the corresponding reduced = 1.102 
(2774.783/2540). While the fit naturally tend towards the afore- 
mentioned two-absorption-l-dust model, we consider additional ab- 
sorbers, driven by our own theory considerations based on thermo- 
dynamic stability arguments (see Section ISTl also Section |4.3. It 
and the de tection of the Un resolved Transition Array in the previ- 
ous XMM dSako et alj|200lh study of IRAS 13349 (Section|4321l. 
While these additional absorbers aid in fitting a few additional weak 
lines, they do not add significantly to the overall improvement in 
fit statistics (AC = 2842.34/2540 = 1.133, = 2741.884/2540 = 
1.093 for the addition of WA-3 and WA-4). - See Tableland Fig- 
ure [13] for details on the four-absorber-l-dust model. We also test 
for the presence of additional absorbing material from within the 
Milky Way and find no strong evidence for its presence along the 
IRAS 13349 hne-of-sight. 

On the intrinsic IRAS 13349 absorption, Figure[T4]shows the 
observed and predicted column densities (A'j ) of the various ions 
we have found in the UV and X-ray spectra of IRAS 13349-1-2438, 
based on both the plasma diagnostic approach of Sections |3] and 
the XSTAR photoionization fitting here. As demonstrated in the fig- 
ure, the series fitting and XSTAR modeling agree with each other 
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Figure 13. The best fit to the Chandra HETGS spectrum of IRAS 13349 (in the observed frame). Numbers associated with the ion labels correspond to 
contributions from one or multiple absorbers as described in Section l43] and Table[8] Note that while this figure is shown with adaptive binning (15 cts bin~^) 
to facilitate easier viewing, our fitting results are based on constant binning at the HETGS resolution. Only predicted lines with equivalent width W\ > 
0.005 mA are labeled. 
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Table 8. Best fit parameters based on fits to the Chandra HETGS data. 
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on the predicted column densities for many of the detected ions. 
For instance, it can be seen that for O VII (Figure \T4\ left panel 
A), the column density predicted from fitting the line series (hori- 
zontal shaded regions) matches with the XSTAR predicted Nj . In 
addition, the figure is instructive in revealing the extent to which 
each absorber is responsible for the different ionic columns. Tak- 
ing the example of O VII again, the dominant optical depth contri- 
bution comes from WA-1 and the Nj due to WA-2 falls short by 
L91 dex. On the other hand, for H-like Ne X, Mg XII and Si XIII 
(panels E-G) the higher ionisation phase of WA-2 satisfies the ob- 
served optical depth due to the respective ion. In the following, we 
consider the impact of adding additional ionized absorbers. 



4.3.1 A third warm absorber? - WA-3 

It is interesting to note that the TVj for He-like Ne IX and Mg XI 
do not fall in the range predicted by either WA-1 or WA-2. This 
makes sense according to panels E and F of Figure [T2I which indi- 
cate that these ions peak at a temperature consistent with an ion- 
isation 2.0 < log^ < 3.0. Indeed, theory considerations based 
on thermodynamic stability arguments (see i]5.1l for details) point 
to more than two discrete zones of absorbers when a 3rd ab- 
sorber (WA-3: logC = 2.98lj;:|;2, Nh = 2.30t;; j2 X iq2i^ ^qjj. 
fined within the ionisation range dictated by our theory consider- 
ations ( i]5.1t is added to the fits, the aforementioned ions are bet- 
ter fit, although changes in fit statistics with the addition of this 
absorber is minimal (AC = 1.143 (2872.129/2540) and reduced 
x2 = L098 (2759.126/2540)). -See Figure[T3]for ion contribu- 



tions from WA-3. We will return to WA-3 in Section ISTI in the 
context of theoretical considerations based on the thermodynamic 
stability of the absorbers in IRAS 13349. In short, the present X-ray 
data are such that the S/N is not large enough for it to drive the fit- 
ting results toward greater than a two-zone absorber, although when 
more ionisation zones are included, the observed column densities 
for H- and He-like Ne and Mg are better matched to that derived 
based on our plasma diagnostic approach of Section |3T| 



4.3.2 UTA Contributions 

As stated, our best model fit to the IRAS 13349 required a dust 
component for the X-ray fitting. For thoroughness, we investigate 
whether additional (gas phase) absorber components can take its 
place. In particular, L- to M-shell photoexcitations in lower charge 
states of Fe I to Fe XVI give rise to absoiption at rest wavelengths 
14 - 17.5 A (~ 0.7 — 0.89 keV), ~ the same spectral region where 
the dust absorption is required. Th ese features, dubb ed t he Unre- 
solved Transition Array or UTA (see Behar et al.l200ll an d lGul2007l 
for theory), was initially detected by|; 



Sako et al. (12001 ) based on an 



XMM RGS study of IRAS 13349. To test for the prevalence and 
strength of these lines based on our method of analysis (i.e. one 
which tie line strengths to the observed UV-Xray ionising SED), 
we explore fits which include additional absorber components that 
may account for them. (As a reminder, note that we began with 
an initial 9-absorber fit, covering the range —4 < log^ < 4, that 
would have accounted for all UTA lines, if significant.) 

We begin by checking the range of ionization parameters 
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Figure 14. The column densities of the vaiious X-ray (panels A-G on the left and N on the right) and UV (panels H-J on die right) found in the high resolution 
spectra of IRAS 13349. Panels K-M, on the right, give the column densities for some of the intermediate theory predicted ions, barely detected in the X-rays. 
The horizontal hashed regions in the left panels A to G, con'espond to derived A'j based on fits to the line series (see Table[3]and Section [3Tt , but weighted 
according to the amount that the X-ray warni absorbers 1 to 3 contributes to it. Respectively, from left to right, green, blue and red points con'espond to the 
best fit A'j as determined by XSTAR fits for absorbers WA-1, WA-3, WA-2. It can be seen that OVII (panel A) is associated with WA-1, whereas ions in 
B-D and E-G panels correspond to WA-3 and WA-2 respectively. In the right panels I and J, the horizontal shaded strips represent the ionic column densities 
A'j derived from the UV observations, assuming a Doppler parameter of 6 = 100 kms~^ (Table|5j. The right panels H-N have a different log^ range, in 
order to accommodate WA-4 and WA-5 (marked with vertical dotted magenta lines). 



where the ion fractions of Fe I to Fe XVI peak (middle panel 
of Figure [TS). As can be seen in the figure, WA-(l-3) primarily 
account for the range of log^ > 1.49, which while fitting for 
UTA ions Fe VII to Fe XVI, do not fully account for the Fe II 
to Fe VI transitions that populate the spectral region between 17- 
17.5 A (0.71-0.73 keV rest; observed ~ 18.8-19.4 A=0.64-0.66 
keV). FigurefTSlshows that the "missing" low-ionization UTA lines 
should fall between —3.0 < log C < 1-0. As such, we include 
two additional warm absorber components: WA-4 (forced within 
-1.0 < log^ < 1.0) and WA-5 (forced within -3.0 < log^ < 
— 1.0) to account for all "missing" UTA producing ions. The re- 
sultant fits (WA-4; log^ ~ 0.99 with logNn ~ -1.78, and 
WA-5: log^ ~ -1.02; with logNn ~ —1.86) primarily affect 
O IV to O VI with the strongest contribution to O V (24.8 A ob- 
served) and O IV (25.2 A observed) absorption, while the fit statis- 
tics (AC = 1.138; xt ~ 1-099) are still poorer than our best 
fit 4WA-l-dust fit, although marginally. (See Section 14.3.31 for de- 
tails on fitting for dust in X-ray spectra.) Indeed, as can be seen in 
Figures[T4land[T5] WA-4 contributes noticeably to the column den- 



sity for these lower ionization oxygen (K-M) and UV ions (H-J), 
while WA-5 give A'j values which are higher than those constrained 
by our HST measurements of Si IV and C IV. Nevertheless, this 
argues strongly for the UV and low-^ X-ray warm absorbers 3-5 
having similar origins, as also borne out based on kinematics (see 
Section |4. 3. 4t . (For subsequent fits, we exclude WA-5 because it 
does not obviously contribute to the fits in a positive way.) It can be 
seen from Figure [T6] (blue) that the UTA alone cannot account for 
all the absorption in this spectral region. 

4.3.3 Direct X-ray detections of iron dust 

As just described, the UTA, while present, is unable to account 
for all absorption in the 17-17.5 A (0.71-0.73 keV rest; ob- 
served ~ 18 .8-19.4 A = 0.64— 0.66) spectral region. It was initially 
proposed by iLee et al" I ( I2OOII) , based on a Chandra HETGS study 
of the Seyfert galaxy MCG -6-30-15, that dust can be directly de- 
tected in high resolution spectra of X-ray bright objects. For iron- 
based dust, these features appear in the form of Fe L (III, II, I) 
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Figure 15. The distribution of ion fractions for the Fe I-Fe XVI ions (mid- 
dle panel), responsible for the UTA, as a function of the ionization param- 
eter For comparison, we also plot the ion fraction distributions of some of 
the relevant, prominent UV (top panel) and X-ray ions (O IV-VII, bottom 
panel). Vertical dotted lines mark the best fit ^ values of the vaiious warm 
absorber components (see Section l4.3.2l for details). 



photoelectric edges between (lab) ~ 14.7 - 17.5 A (0.7-0.84 keV; 
overlapping with the UTA spectral region). (Fe-based dust in astro- 
phy sical environments can also be measured at ~ 7 keV Fe K - 
see iLee & Raveil 12001 but this is not relevant to the Chandra 
HETGS capabilities.) As can be seen in Figure[T6](blue), the UTA 
is unable to account for much of the absorption between 19.3- 
19.5 A (observed; ~ 0.63 - 0.64 keV), consistent with the Fe L III 
edge (redshifted to 2/^13349) - in the context of IRAS 13349, this 
argues strongly for direct detections of iron-based dust in the AGN 
environment. 

To investigate the direct detection of iron-based dust in 
IRAS 13349, we incorporate cross sections associated with con- 
densed matter forms of iron into our 4-absorber fit to determine 
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Figure 16. The Chandra MEG observed spectrum zoomed into the UTA and 
iron dust spectral region, at the redshift of IRAS 13349 (black). It is cleai' 
that ionised gas alone (5-warm absorbers including the UTA) is unable to 
explain all the absorption in this region (blue; see Section l4.3.2l for detailed 
discussion on UTA), and additional absorption with origins in iron-dust in 
IRAS 13349 is also needed (red; Section l433t . 



an ionic column for iron in grains to be A^Fcdu 



(1.2 ± 



0.21) X 10^^ cm~^ (See Lee et all l2009l for details on X-ray 
methodology for determining the quantity and comp osition of in- 
terstel lar dust.) Using the ISM abundance values of I Wilms et al.l 
ll2000h . this translates to an equivalent Hydrogen column A^h ~ 
(4.5 ± 1.30) X 10^^ cm"2. While the X-ray data are not of suffi- 
cient S/N to distinguish the exact iron composition (e.g. pure Fe vs. 
FeO vs. Fe2 03 vs. Fe2 04), a conservative estimate of an observ- 
able transmission signal, requires that the iron-based grains have 
thickness t ~ 0.1 — 0.8 pim, based on the formalism for transmis- 
sion T = where p is the density of the specific compound 
as defined from "The Handbook of Chemistry and Physics"; the 
values for the attenuation length were obtained from theFfcXRO at 
LBL, and[£|NIST. Furthermore, the fits strongly suggest that at least 
90% of the Fe is locked up in grains, with a minor contribution from 



gas phase iron, A^Pe- 



2.1 X 10^^ cm" 



4.3.4 The absorber kinematics 

Of interest, the outflow velocities for the X-ray absorbers we detect 
are best matched to the higher velocity (~ 950 km s^^ ) UV out- 
flowing component. -See e.g. Figure[T7]comparison of the HETGS 
view of O VII versus STIS view of N V. Figure [14] shows that the 
observed column densities for the UV ions C IV and N V (panels I 
and J) are consistent, within eiTors, with the XSTAR predicted col- 
umn densities for the low ^ phase of the warm absorber, determined 
from fitting to the Chandra spectra (OVII, panel N). Our analysis 
shows that the low-^ (WA-1 and WA-4) can have linked UV and 
X-ray absorption from C IV to O VII and some O VIII, although the 
latter is primarily associated with the higher-^ WA-2 and WA-3. 

Furthermore, it would appear that the lower ionisation phases 
are associated with a higher outflowing velocity. Our lower limits 
on the derived values for Ubiuo are roughly consistent with the up- 
per limit findings by Sako et al., (2001; «biue = 420l kms"^ ) 



° http://www-cxro.lbl.gov/opticaLconstants/ 

^ http://physics.nist.gov/PhysRetT)ata/FFast/html/form.html 



© 0000 RAS, MNRAS 000, 000-000 



22 



Lee et al. 



o 

CO 
m 

'a 

6 

u 
o 
c 

CO 
X 

3 



VII zi.6o;i,' is'.eB?,' iv.ves! 17.396, i^'aoo 










MM 


.A 

MM 



N V 1248.80, 1238.82 




i I . . I 



-4000 -2000 2000 

Velocity 



4000 



Figure 17. A velocity profile comparison of the X-ray and UV absorbers, 
based on representative ions show that the UV and low-^ X-ray warm ab- 
sorbers likely have a common origin. ("TOP): velocity spectrum combined 
from the five strongest (X-ray) resonance lines of O VII based on Chandra 
MEG X-ray spectra, and (BOTTOM) velocity spectrum combined from two 
strong N V absorption lines seen in HST STIS spectra. Details for how these 
are generated can be found in Section lTTI 



and Holczer et al, (2007; ?;biuc ~ 300 ± 50 kms ^ ) ;/ we 
account f or the difference in choice of redshift value, i.e. z — 
0.10764 jKim etalj[T99l) used by previous authors versus the 
z — 0.108530 value derived in this paper based on higher 
resolution HEX data. We do not find any evidence for the ~ 
+20 kms^^ outflow reported by Sako et al. however, and our fit- 
ted turbulent velocity widths are slightly lower than the Wturb 
640 km values reported by both Sako et al., and Holczer et al. 
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Figure 18. Stability curves illustrating differences in the thermodynamic 
behaviour of the absorbing gas as illuminated by the full UV-X-ray SED 
(solid red curve conesponding to the SED in solid red shown in Figure [Tol 
bottom panel), versus an X-ray only ionising SED (black dashed curve cor- 
responding to same representation of ionisi ng SED in Fig ure \iQ\ bottom 
panel). The unstable regions identified by the lHolczer et alj ^2007) analysis 
of this Chandra X-ray data based on an X-ray only ionising SED is high- 
lighted in bold black. Also in bold on the red curve is the aforementioned 
unstable region translated over the same AT for the stability curve drawn 
based on a UV (i.e. accretion disc) inclusive SED. The filled circles on the 
red curve are the points con'esponding to WA-1, WA-2, WA-3, WA-4 and 
WA-5 as predicted by the XSTAR fits to the Chandra HETGS data and the 
HST STIS data and/or theoretical considerations, as discussed in Section 
14.31 As can be seen, no unstable regions are identified by our analysis. The 
dotted blue lines connect the points with the same ^ values in the two curves, 
thus showing that the temperature for all the ^ values are higher for the gas 
ionised by the X-ray only SED, except for WA-5, where the temperatures 
are equal. 



5 DISCUSSION 

5.1 Thermodynamics and structure of absorber 

A topic of active debate is whether warm absorbers in AGN are in 
a continuous medium or are an ensemble of discrete clumpy me- 
dia which are in different thermodynamic phases in near pressure 
equilibrium with each other. The answer to this question has in- 
teresting consequences for the geometry of the AGN environment, 
as a whole. Here, we use stability curves to demonstrate that the 
choice of ionising SED (e.g. X-ray versus the full source SED) for 
IRAS 13349 is very relevant in the conclusion for one scenario over 
the other. 

Stability curves are thermodynamic phase diagrams of tem- 
perature (logT) versus pressure (log(5/r), where ^ = ionisation 
parameter. It is an effective theoretical tool often used to discuss 
the thermodynamics of p hotoionized gas as s ociated with the X-ray 



and UV absorbers (e.e. iKrolik et al.l 19811 


Krolik(&Krissll200ll: 


iRevnolds & Fabianlll99a Chakravortv et al. 


20091 and the refer- 



ences therein). By definition, an isobaric perturbation of a system 
in equilibrium is represented by a small vertical displacement from 
the stability curve; such perturbations leave ^/T ^ L/pR^ con- 
stant, which for constant L/R^ leave the pressure unchanged. Any 



'system' located on the part of the curve with positive slope is a sta- 
ble thermodynamic system, because a perturbation corresponding 
to an increase in temperature leads to cooling, while a decrease in 
temperature leads to heating of the gas. If the stability curve is char- 
acterised with discrete allowed values of temperature at the same 
pressure, it points to a cloud-like absorber, whereas a wind-like 
scenario will have a continuous distribution of allowed temperature 
and pressure. Before proceeding, it is informative to briefly discuss 
the definition of the ionisation parameter in the context of stability 
curves. In the definition ^ ~ L/nR^ , it is conventional to use the 
luminosity L in the energy range 13.6 eV — 13.6 keV (i.e 1 — 10^ 
Rydberg). The warm absorber properties are, however, determined 
by the photon distribution in the soft X-rays {E > 100 eV) and 
not by photons with energy E « 100 eV. Hence, in the mak- 
ing of warm absorber stability curves, authors often use modifi- 
cations in the definit i on of the ionisation parameter. For example, 
IChelouche & Netzed ( [200l and references therein) use an ionisa- 
tion parameter Ux which considers the ionisi ng flux only between 
540 eV to lOkeV. IChakravortv et al.l ( l2012h . on the other hand, 
simply normalise the stability curves by using the ratio of the total 
luminosity to the luminosity in the X-ray power-law component in 
the SED. However, it is to be noted that these modifications simply 
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shift the entire stability curve along the log(^/r) axis and does not 
change the nature of the kinks or curves in the stability curve. For 
example, the range (A[log(^/T)] or AT) of the stable or unstable 
regions of the stability curves would remain unchanged with such 
modifications of the ionisation parameters. As such, all the qualita- 
tive properties of the stability curves and comparisons between the 
stability curves of different SEDs, discussed subsequently, would 
remain the same even if better definitions of the ionisation parame- 
ter are used. 

In the context of AGN warm absorbers, several authors (e.g . 
iKrolik et alJl98lMKrolik & Krisj200ll ; lGehrels & William j 19931) 
have commented that the gas in the 4.5 < logT < 5 region is 
thermally unstable, since the cooling function A(r) has a nega- 
tive gradient in this temperature range. As such, the ions detected 
on the stable regions of the stability curve bracketing this temper- 
ature range would be attributed to cloud-like absorbers where the 
different phases are in pressure equilibrium with each other. Based 
on the findings for j ust such an ionisation gap in IRAS 13349 by 
lHolczeretalT ( l2007h using data from the Cha ndra HETGS (the 
same data we are re-analysing in this paper) and lSako et al.lJlOOlh 
with the XMM-Newton RGS, cloud-like clumpy absorbers would 
appear to be the c onclusion one wo uld draw for this quasar. (We 
note however that ISako et al. I I2OOII acknowledge that their stabil- 
ity curves do not show thermodynamic ally unstable regions for 
IRAS 13349.) In a comparison with the iHolczer et al.l | |2007|) re- 
sults, the fact that their absorption measure distribution (AMD) 
methodology, despite assuming a continuous distribution of ions 
at the onset, concludes for thermodynamically unstable regions, 
further strenghtens the case that independent of fitting technique 
(AMD vs. photoionisation modelling as that employed here e.g.), it 
is the ionising SED that is the critical factor. 

.Revnold s & FabiarJ l ll995h and lChakravortv et al.l ( l2012h dis- 
cuss and demonstrate that a moderate to strong soft excess com- 
ponent makes the s tability curve more stabl e at 10^ K (see also 
lKroliketaLlll98ll) . IChakravortv et all ( l2012h further demonstrate 
in general for AGN that the inclusion of the accretion disc spec- 
tral component in the EUV range in addition to a significantly 
strong soft excess may completely remove the apparent thermal 
instability in the aforementioned temperature range. In this paper, 
we revisit the conclusions for the thermodynamic stability con- 
ditions for the IRAS 13349 warm absorber gas in the context of 
different SEDs; one that excludes the disc components (aforemen- 
tioned X-ray studies of IRAS 13349), and one which includes them 
(this paper). For our simulations of the photoionized gas thermo- 
dynamic properties, we used XSTAR to model a spherical shell of 
solar metallicity gas having A^h = 5 x 10^^ cm^^ and particle den- 
sity riH = 10^ cm"'*, illuminated by the aforementioned SEDs. It 
should be noted that the stability curves are insensitive to the choice 
of 10^ < n < lO^^cm"^ for the IRAS 13349 ionising spectra. It 
can be seen from Figure [TS] that the choice of ionising spectrum 
strongly influences the conclusions for gas thermodynamic stabil- 
ity. For the two curves shown, we use thicker lines to highlight the 
temperature range conventionally established to be unstable (neg- 
ative slope) for photoionized gas in AGN environments. It can be 
seen that an X-ray only (disc-free) ionising SED leads to the pre- 
diction (black dashed curve) for a narrow range of ^ coincident 
with the aforementioned temperature range over which the warm 
absorber is mildly unstable. Indeed, it is exactly the finding for ions 
on eith er side of this unstabl e phase seen in the black dashed curve 
that led lHolczer et alj ( |2007|) to their conclusions for IRAS 13349. 
In contrast, when we include the UV as part of the ionising input, 
no such instability is seen, thereby pointing to a flow that is more 



wind-like and continuous in IRAS 13349. The effect of the excess 
UV photons due to the addition of the accretion disk component, 
facilitates the formation of OVII and other ions with similar ion- 
ization potential (see panel D of Figurell2t. OVII acts as a cooling 
agent and hence lowers the temperature of the gas, which results 
in thermodynamic stability in the region of WA- 1 . Moreover, the 
accretion disk component also helps to lower the Compton tem- 
perature of the gas. Hence the temperature of the curve in the re- 
gion of WA-3 is lowered, rendering this part of the curve stable, 
as well. Thus, while our analysis reveals similar ionisation states 
found from previous studies, the additional consideration of the 
UV-influence drives us to a conclusion that the ionised absorber 
in IRAS 13349 should be more continuous and wind-like than dis- 
crete clouds. 

Taking the thermodynamic instability-free red curve in Fig- 
ure [TS] it is interesting to further consider the conditions driving 
the fit results. Specifically, as discussed in Section |43] the compo- 
nents WA- 1 and WA-2 are statistically 'sufficient' to fit the Chan- 
dra HETGS data. However there are no thermodynamically unsta- 
ble phases between these two points to separate them — as seen in 
the figure (red curve). In other words, there is a range of ^ which is 
theoretically allowed by thermodynamic calculations and yet is not 
required while fitting the data. While one explanation for this may 
be the insufficient S/N of the available data, there may be additional 
more physically motivated answers. Specifically, (1) there may be 
an absence of absorbing material of sufficient column density at 
certain intermediate ionisations (which may be related to interme- 
diate distances from the black hole) to imprint significantly on our 
spectrum. Indeed, for IRAS 13349, when we 'force fit' a third com- 
ponent WA-3 to the data, its physical conditions (^, Nh and Vbiue) 
come out to be intermediate between that of WA- 1 and WA-2. Fur- 
ther, (2) when we include additional theoretically motivated com- 
ponent WA-4 and WA-5 to account for the entire UTA, we indeed 
populate the entire stability curve suggesting presence of material 
in a continuous distribution of ionisation states The stability condi- 
tions of a super-solar warm absorber would be distinctly d ifferent 
from that o f a gas with Solar abundances (see Chakravo rtv et al.l 
l2009ll2012h . Since the outflowing gas from the nuclear regions of 
AGN may indeed have super-solar metallicities (e.g., lArav et aU 
we explore this idea further for IRAS 13349 in the next sec- 
tion. 

5.1.1 Super-solar metallicities 

IChakravortv et al.l ( l2009l l2012h show that the stability curves for 
a super-solar absorber is more likely to show discretely separated 
multi-phase structure. In Figure [19] we show the stability curves 
for different super-solar metallicites of the absorber. Independent 
of whether only iron, or all metals are varied, it can be seen that the 
higher metallicity gas starts exhibiting thermodynamic instability 
at logT > 5.2, and the high ionisation phase WA-2 (log^ ~ 3.4; 
top-most, green circle) falls in a part of the curve after it recov- 
ers from this instability. However, there is no distinct unstable re- 
gion between the low WA-1 (log^ ~ 1.7) and intermediate WA-3 
(log C ~ 3) ionisation phases which would mean a continuous dis- 
tribution of phases in this range of ionisation parameters. WA-3 
falls exactly at the turnover where the high metallicity curves be- 
come unstable at logT > 5.2, suggesting that while WA-2 and 
WA-3 are thermodynamically separated, they are nevertheless in 
pressure equilibrium. It is of further interest to note that the WA- 
2 and WA-3 values derived from XSTAR fits to the Chandra data 
fall on opposite ends of the unstable regions in Figure [T9l-This 
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Figure 19. Stability cui'ves corresponding to super-Solar gas of 3Z0(blue), 5Z0(pink), 7ZQ(black), and 9ZQ(green), where (LEFT) Zq refers to all 
Z-elements, and (RIGHT) Zq refers only to iron, i.e. only iron abundance are altered here. For compaiison, in both panels, solar abundance stability curves 
are shown in solid red, and plotted points correspond to XSTAR predicted ionisation phases for the warm absorber 




Broad Line Region 



5.2 The IRAS 13349 dusty warm absorber 

Given the iRevnoldsl l ll997h finding back in the ASCA era that 
~ 20% of the warm absorber Seyfert Is have significant optical 
reddening, it is not surprising that there is now a remarkable con- 
vergence of evidence for dusty ionised winds from observations 
with more powerful satellites. T hat IRAS 13349 harbou r s dust has 
long been establish ed with IR iBeichma n et al 19861: iBarvainisI 



1 19871). UV/optical jWillsetidJ I199Z ; Hines et alj l200lh and X 
ray iBrandt et all 1 19961 1 19971) studies. Using the superb multi 



wavelength data set presented here, in combination with knowl- 
edge based on prior work, we present a geometrical model for the 
IRAS 13349 line-of-sight (LOS), and location of dust, based on 
dust modelling of the absorbed optical-UV spectra (Section[5]2j2} , 
and direct IR and X-ray measurements (Section|5]23} of dust prop- 
erties. 



Figure 20. Inferred geometry of IRAS 13349. See '^^tor details. 



may not be coincidental, since gas in the unstable regions would 
preferentially migrate towards stable region, by cooling or heating. 
Ultimately, if there are discrete (rather than continuous) absorbers 
in IRAS 13349, at least some components of the gas must have 
super-solar metallicities. 

We would like to caution that our discussion of super-solar 
metallicities is not quantitatively absolute and should only be con- 
sidered as a qualitative impression of our expectations. A more rig- 
orous study would involve fitting the data with super-solar popula- 
tion files in XSTAR and drawing conclusions based on the best-fit 
metallicity of the absorber. Given the S/N of our data, however, 
such detailed fits are not warranted in this paper 



5.2.1 A model for the IRAS 13349 dust 

As we showed in i|4.1| the direct UV and optical flux from 
IRAS 13349 is substantially lower than would be expected for a 
typical QSO. Additionally, the high polarisation rising to shorter 
wavelength s llHines et alj|200lh is indicative of scattering by dust. 
iHines et aljj200lh attempted to model the UV/optical flux and po- 
larisation of IRAS 13349 as a combination of a direct view of the 
nucleus reddened by an intervening screen, plus dust-scattered light 
from biconical regions filling the poles above the obscuring toius. 
Using dust with SMC-like properties, they were able to obtain a 
reasonable match to the shape of the polarised flux spectrum, but 
the predicted total flux spectrum did not match the observations — 
the optical flux was much too high, and the reddening caused the 
predicted UV flux to drop far below the observed intensity. When 
one also considers the shorter-wavelength data from STIS pre- 
sented in this paper, the fact that the UV flux is much higher than 
predicted by their model suggests that there could be additional 
scattering that is colour neu tral, e.g., like the Thomso n mirror of 
hot electrons in NGC 1068 jAntonucci & Millei|[l985h . This idea 



© 0000 RAS, MNRAS 000, 000-000 



IRAS 13349+2438 25 



is reminiscent of the original suggestion bv lWills et al.l ( [l993) tliat 
tiie polarised light in IRAS 13349 was due to Thomson scattering, 
although that idea proved inconsistent with the waveleng th depen- 
dence into the near-UV observed by iHines et al] ( 1200 ll) . For our 
proposed geometry here, the Thomson-scattered component would 
play a minor role in the polarisation, and only account for the flux 
at the shortest wavelengths. 

To test this idea, we developed a simple heuristic model as 
illustrated in Fig. |20] As shown in the figure, our direct line of 
sight (the dash-dot line) to the active nucleus passes through the 
upper atmosphere of the obscuring torus. The atmosphere consists 
of a mixture of highly ionised gas and dust that has not yet been 
evaporated by the incident radiation. This line of sight is heavily 
reddened in the UV/optical band, or even completely blocked; in 
the X-ray, since the gas is highly ionised, there is litt le absorption 
by lig hter elements, and the active nucleus is visible jSrandt et al.l 
Il997h . A second line of sight (the solid line and the dashed line) 
passes through less dense portions of the atmosphere that are more 
transparent at UV/optical wavelengths, but this line of sight only 
gives a view of the active nucleus as reflected by dust from the 
far wall of the obscuring torus. Finally, a third line of sight passes 
above the atmosphere of the torus, and it has a view of the nucleus 
that is Thomson scattered from hot electrons in the thermal wind 
escaping along the polar cavities of the torus. 

Assuming the active nucleus plus BLR has a spectrum identi- 
cal to the composite QSO spectrum (defined as fc), this model can 
be analytically expressed as: 

i^optuv(A) = [Fdircct(A) + F,cat(A) + FTb{X)]xlO-°-^^'''^^\5) 

where Av is the optical absorption within our Galaxy and 
g{\) — Ax /Ay is the ex tinction law in the Milky Way, as 
defined in lFitzpatric3 ( Il999l) . The other components within the 
IRAS 13349 galaxy are: 

• Absorbed direct flux in the torus, torus atmosphere, and 

galaxy: We define the LOS flux passing directly through the torus 
by: 

fdircct(A) =rd/c(A) (6) 

where Fd is the absolute normalisation of the composite spectrum, 
Avtor ™d (7tor(A) = ^Ator/^Vtor t^e optlcal absorption and 
the wavelength-dependent extinction law within the torus, Ay^tn, 
and (7atm(A) = ^Aatm/^Vatm ^6 the optical absorption and the 
wavelength-dependent extinction law within the torus atmosphere, 
and ^Vgai and (7gai(A) — ^A^^i/^Vg^i are the optical absorption 
and the wavelength-dependent extinction law within the galaxy, 
respectively. 

• Absorbed dust-scattered UV/optical flux in the torus atmo- 
sphere and diffuse medium: 



• Absorbed Thomson-scattered UV/optical flux in the diffuse 
medium: 



Fecat(A) =7.catrd/c(A) ( 1+^ 



(7) 



xlO 



-0-4 (-^V^tm 9atm(A) + Av„, 9gal(A)) 



where (^1 + x") stems from the wavelength dependence of the 
scattering (a = 4 is the Rayleigh limit), Ac is the critical wave- 
length of the scattering, characteristics of the typical dust grain 
size, and 7scat is the fraction of direct light which is dust-scattered. 



-Prh =7Thrd/c(A) X 10 



-0.4 A^, 



(8) 



where 7Th is the fraction of direct light which is Thomson scat- 
tered. 



5.2.2 The wavelength-dependent extinction within IRAS 13349 

With this model in mind, we consider, through fitting, which 
combination of dust absorption and scattering best-describes 
the IRAS 13349 obseri'ed optical-UV spectrum. We used the 
Levenberg-Marquardt minimis ation to determin e goodness of fit, 
based on the IDL code MPFIT ( lMarkward3l2009h . However, we do 
not have any information about the extinction within IRAS 13349, 
neither in the galaxy nor in the torus and torus atmosphere. 
We therefore considered several p ermutations of ex isting laws 
previously published for the MW (lFitzpatricklll999h . the aver- 
age SMC and LMC dOordon et alj|2003h . and starburst galaxies 
l iCalzetti et al.ll2000h . The worst fits were obtained with the MW 
and/or LMC extinction, due to the presence of the 2175 A bump 
in these laws. This points towards IRAS 13349 lacking the un- 
known dust population responsible for this 2175 A absorption fea- 
ture in the MW and the LMC, in good agreement with several stud- 
ies showing that the dust content in AGNs is SMC-like (see e.g. 
ICrenshaw et aflEoOll : lHopkinsll2004l and references therein). As 
expected, the use of the SMC and/or the Calzetti laws gives better 
results, but the fits are still poor, with reduced ^ 3. 

We therefore build our own SMC-like extinction laws for the 
IRAS 13349 torus, torus atmosphere, and galaxy as: 



5tor(A) = 



in the torus. 



5atm(A) 



/ 5448 A 

/ 5448 A 
in the torus atmosphere. 



.9gal(A) = 



/5448 A\'^'=''' 



(9) 



(10) 



(11) 



in the galaxy. These laws are built so that they are equal to I in 
the V band, and we assume the dust population is different in 
the torus, the torus atmosphere and the rest of the galaxy. In our 
fit, the spectral indices /3tor, /3atm, and /3gai are free parameters 
along with the optical extinction within the torus Avt„, the 
atmosphere Av^tm ' ™d the galaxy ^Vg^i , the fraction of dust- and 
Thomson-scattered light 7scat , 7Th, and the critical dust-scattering 
wavelength Ac. The absolute normalisation of the composite 
spectrum Fd is fixed to the level discussed in Section 14.11 the 
dust-scattering spectral index a is fixed to 4 (Rayleigh limit), and 
the optical extinction within the MW is fixed to Avq^, = 0.04 ala 
the lWills et ai] ( ll992h derived value for IRAS 13349." 

• Best-flt dust model: In our fitting, we are unable to constrain 
the fraction of dust-scattered light 7scat as it systematically pegs 
to 1. We can interpret this as due to (1) an overestimation, by the 
composite spectrum, of the nucleus continuum in the optical or, 
(2) the presence of another contributing component in the opti- 
cal domain. We therefore fix 7scat = 0.2, a conservative value 
chosen to reflect the absolute normalisation of the absorbed UV 



© 0000 RAS, MNRAS 000, 000-000 



26 Lee et al. 



Frequency (10 Hz) 
20 15 



10 



Observed spectrum 

Total light ■ 
Dust-scattered light ■ 
Thomson-scattered light 
Polarized emission ■ 




1 2 3 4 5 6 7 8 

Observed wavelength (1 0'^ A) 

Figure 21. Best-fit model for the UV/Optical spectrum of IRAS 1 3349 (ob- 
served spectrum in grey). The full modelled spectrum is the red hne. The 
Thomson-scattered component is the green line, and the dust-scattered com- 
ponent is the blue line. We also include the polarised light (magenta) as 
measured in Hines et al. (2001). Our fit results point to the direct light hav- 
ing marginal contribution, pointing to almost complete obscuration by the 
torus in the UV and optical. 



spectrum needed to bring IRAS 13349 in line with unabsorbed 
QSO behaviour (see Section 143] for details). Likewise, it appears 
difficult to constrain any torus parameters, since all fits tends to- 
wards removing the direct light. This strongly suggests that the 
direct light is completely blocked by the torus, and so we fit us- 
ing only the dust- and Thomson-scattered light. Our best fit (Fig- 
ure |2T] and Table |9} points towards a steep extinction law within 
the torus atmosphere {Paxm ^ 3.01), while the extinction within 
the galaxy is flatter, with /3gai ~2.11, but steeper than that of the 
SMC. The typical grain size derived from the dust-scattering func- 
tion a ~ Ac/27r ~ 1.21 ^m, is consistent with large grains within 
the torus atmosphere absorbing the dust-scattered light in the far- 
UV, where the emission is dominated by Thomson-scattering, with 
7Th ~ 2.4%. We have checked that the lower normalisation for 
the UV peak (~factor-of-two compared with the .Elvis 1994 value 
we use ; see Section l4~Tl for details) predicted by the lHopkins et al.l 
( |2007|) SED give roughly similar dust fit results. 

We attempted to, but ultimately could not include polarisation 
in any reasonable way in these fits to gamer believable measure- 
ments, although the conclus ions based on our p olarisation-free fits 
are consistent with those of iHines et alj ( 1200 ll) based on polarisa- 
tion studies. Specifically, Hines et al. showed that polarisation is 
wavelength-dependent, therefore pointing to dust-scattering dom- 
inating over wavelength- /^dependent Thomson-scattering. Simi- 
larly, our fits point strongly to dust-scattered light dominating the 
2000 and 8000 A spectral region. 



5.2.3 Direct X-ray and IR measurements of the IRAS 13349 dust 

Using additional knowledge gained from direct X-ray and IR 
measurements of the IRAS 13349 dust properties as enabled by 
the superb spectral resolution of the Chandra HETGS and Spitzer 
IRS, we further consider the location of the different composition 
grains we detect. 

• Chandra X-ray spectroscopic detection of iron grains: 

As described in Section 14.3.31 we directly detect iron-based 
dust in the X-ray corresponding to Nup^.Q^^i^^ = (4.5 ± 1.30) x 



10^^ cm ^ and grain size ~ 0.1 — 0.8 fim .According to the 
radiative transfer calculations of iKama etaL I ( l2009l) . iron subli- 
mates between ~ 1000-2700K as a strong function of densities 
(i.e. pressures). The protoplanetary disk conditions that they con- 
sider has typical vapour densities of ~ 10^^'^-10~^''g cm"'* which 
translates to particle density ~ 10^*^-10^^ atoms cm~^, assuming 
55 grams moleFe"^, not dissimilar to the particle densities that 
one would expect for the inner AGN clouds. More importantly, at 
these densities, iron sublimates at 1400K, a higher temperature than 
silicates at 1 100-1300K. As such, in the context of our geometrical 
model for the IRAS 13349 viewing geometry (Figure [20). we can 
place the iron in the upper, less optically thick layers of the torus 
atmosphere, with the silicates more towards the interior Further- 
more, in adopting the formalism 
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Vl045ergs-V V Tsub j 

derived by iNenkova et al. I ( l2008h , we find ^ 2.7 — 5.2 pc for 
the inner radius of cloud distributions in a clumpy torus for sub- 
limation temperatures T~ 1100 — 1400K and Lboi ~ 3 x 
10 ^^ erg s~^ (see Table 7) . The Keck Interferometer measurements 
of iKishimoto et alj ^20091) point to dust at ^ 0.92 pc, which sug- 
gests either dust T ~ 2000K or larger dust grains (e.g. 1.2/^m 
as derived based on spectral fitting in Section 15.2.1) - see also 
l lSitkoetalJl9"93l e.g.). 

• Spitzer IRS spectroscopic detection of silicates and PAHs 

As discussed in Section 13.41 the Spitzer IRS spectra show 
weak 7.7 and 11.3 /xm PAH, as well as broad 10 pan and 
18 pan silicate emission. The PAHs are associated with the outer 
starbursting part of the galaxy and therefore play no role in the dust 
absorption and scattering effect considered previously. Even if the 
obscuration is a supernova heated starburst disk (see e.g. theory in 
Fabian et al. Ill998l : IWada & Normanll2002l : iThompson et"ai]|2005l: 



IBallantvnell2008h . in the IRAS 13349 host galaxy, any PAHs will 
have to be far enough away from the strong radiation to have little 
effect. Else, for the detected silicates, given that they are in emis- 
sion, just like in other normal quasars and Seyfert Is, it is quite 
clear that we are not looking through a lot of dust, but rather are at 
a high-enough elevation angle above the torus to clear the densest 
gas, but perhaps just grazing the surface and looking through its 
atmosphere, which is what is causing the UV and X-ray absorption 
lines and the extinction of the UV continuum. 

Furthermore, the lO^m peak is shifted to the red, indicating 
either a depletion of small dust grains or more crystalline sili- 
cate dust (e.g. enstatite) rather than amorphous olivine. The opac- 
ity at the 1.844 keV (6.724 A, rest) Si K edge is too low for 
us to verify the type of silicates with X-ray data. However, the 
fact that the dust we do detect in the X-rays is primarily locked 
up in iron , argues strongly for the silicat es to be Mg-rich (e.g. 
forsterite). iMarkwick-Kemper et al.l ( l2007h argue for forsterite in 
the PG 2112-1-059 BAL QSO wind based on Spitzer studies; one 
might imagine a similar dusty outflow for IRAS 13349, but as 
stated in Section [3.4.2l assessing the level forsterite may contribute 
to the 1 1.2 feature is beyond the scope of this paper 

In further estimating where the silicate dust responsible for 
the broad IR emission may reside, we derive fisub ~ 1-2 pc, us - 
ing the relation 7?sub = 0.2(Lboi/10'"^)^/^ pc ( iLaoili 19981^0031) . 
To arrive at this, we take the mid-IR luminosity of IRAS 13349 
i/L^(25/im) = 3.2 x lO"'' = Lboi/c/i erg s"\ where /c ~ 0.6 
is a representative dust covering fractio n and /j '-^ 0.5 is a repre- 
sentative mid-IR anisotropy for quasars ( lOgle et al .120061 ; [s hi et al.l 
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I2005I) . Based on this, it is likely tiien that the dust responsible 
for the mid-IR continuum is the same dust that absorbs the direct 
optical-UV spectmm, since it accounts for a large fraction of the 
absorbed and reprocessed bolometric luminosity of the quasar. 

5.3 The dusty wind in IRAS 13349: torus or accretion disk 
origin? 

Given the excellent qualitative agreement between our proposed 
model for the dust contributions (Section |5.2.2t and the observed 
spectrum, our picture in Figure |20]must be close to reality. One ad- 
ditional piece of evidence in suppo rt of our model can be gleaned 
from the existing polarisation data. iHines et all l l200ll) note a dif- 
ference in polarisation angle for the broad emission lines com- 
pared to the continuum. The line centres have their angles rotated 
by 5 deg, but the polarisation stays fairly constant. One can inter- 
pret this rotation in terms of the differing angle between the scat- 
tering region and the sources of the incident radiation, meaning 
that the BLR is displaced from the optical/UV continuum source 
by several degrees as viewed from the scatterer. Also, since the 
emission lines are not depolarised, the angles from across the full 
BLR to the scatterer must be within a fairly naiTow range so that 
their polarisations do not cancel out. So, the size of the BLR as 
seen from the scatterer must be small. If we take the difference 
in polarisation angle as the offset of the BLR from the contin- 
uum source, then, given a distance for the BLR from the nucleus, 
we can deduce the distance of the scatterer from the continuum 
source. Scaling up reverberation mapping results for BLR sizes 
jKaspi et alj|2005l) to our extinction-corrected UV luminosity for 
IRAS 13349 (L1350 = 7.95 x 10'*'' ergs"^ ), the BLR would be 
at a radius of r = 7 x 10^^ cm. The scatterer would then be at 
8 X 10^* cm, or 2.6 pc. For comparison, the dust sublimation ra- 
dius Rsuh ~ 1. 2 pc we calculate in Section |5.2.3l (same value from 
lBarvainij[l987h . and the inner ed ge of the torus is expected to be 
at ~ 1 pc(L/10''^)°-'^ = 2.8 pc iKrolik & Begelmaij|l988h . This 



is remarkably consistent with the optical light being scattered light 
from the far wall of the torus. In adopting the classic picture pre- 
sente d whereb y an obscuring torus is a lways present dAntonuccil 
1 1 9931 : see also iLawrence & Elvis|[T982h our view of IRAS 13349 
then is at high inclination, barely skimming the surface of the toms. 

This picture presents a problem then for scenarios of 
accretion-disk winds (e.g.. lElvisll200Clh . where the wind stream- 
lines favour the equatorial plane, parallel to the surface of the torus. 
Nevertheless, it is of interest to consider further an alternative to 
the torus for the obscuration, that would allow the wind t o have 
an accretion disk origin. Recently, ICzemv & Hrvniewiczl ( l20Ilh 
proposed a model for the creation of dusty winds at > 10007?s 
by drawing parallels between clouds on quasar disk surfaces to 
similar temperature and pressure conditions that cause dust for- 
mation in AGB stars. In the context of the Elvis quasar model, 
the material in this wind would then be the cause of the obscu- 
ration, in lieu of the torus. But can this satisfy our condition for 
a high inclination viewing angle for IRAS 13349? Based o n the 
wind flow lines shown in Fig. 3 of iRisaliti & Elvi j | |20IO|) . the 
scale height at which our line-of-sight intersects the dusty wind 
isz = 2007?s = 0.09 pc at lOOOiJEg ^ for a 10^ Mq, quite close 
to the dis k plane. Alternatively, a w arped disk as the source of ob- 
scuration jLawrence & Elvisll20Iol) . although for IRAS 13349, it 
would have to be contorted in such a way as to satisfy the known 
effects caused by dust scattering, absorption, and polarising. Else, 
if we again consider the putative torus, one might argue that our line 
of sight on the near side of the torus passes under the streamlines 



Parameters 


Best-fit value 


Definition 


^scat 


0.2 (fixed) 


Fraction of dust-scattered light 


Ac (fJ-m) 


7.60 ±0.62 


Dust-scattering critical wavelength 


a 


4 (fixed) 


Dust-scattering index 




0.14 ± 0.02 


Optical extinction in atmosphere 




3.01 ± 0.08 


Spectral index of extinction in atmosphere 


7Th 


0.024 ± 0.003 


Fraction of Thomson-scattered light 


Ay , 


0.05 ±0.01 


Optical extinction in galaxy 


/3gal 


2.11 ±0.11 


Spectral index of extinction in galaxy 


X? (d.o.f) 


1.66 (862) 





Table 9. Best-fit parameters. Uncertainties are given at the 2a confidence. 



of the accretion disk wind, but then the line of sight from the nu- 
cleus to the far wall of the torus and back again must pass directly 
through the wind, and at multiple angles. 

While the accretion disk wind idea is a desirable one, it is ul- 
timately more difficult to reconcile, compared to a torus wind, for 
IRAS 143349. The line widths in our STIS spectmm are similar 
to the narrow, few hundred kms"*^ absorption lines common in 
other nearby AGN jCrenshaw et"aLll2003l: lKrissll2002l: iDutiii et al.l 
' 20071), as expected for the thermal winds of Krolik & KrissI 1 19951 
120011) . Also, the velocity shifts we detect in both UV and X-ray are 
slow compared to th e ultra-fas t outflows expected from the disk 
jTombesi et alj l20ld : see also IChartas et atl |2002| : IPounds et al] 
I2OO3I for the first claims of such outflows in Seyferts). 

If there is no disk wind in this object, then perh aps the condi- 
tions t h at drive it are just not righ t. The disk winds of lMurrav et al.l 
( ll995h : lMun-av & Chiang| ( ll995h require a high ionisation parame- 
ter and and a high column density at the inner edge of the wind, and 
IRAS 13349 does not have a high column density of either X-ray 
or UV-absorbing gas. The manifestation of a disk wind in an AGN 
may not be a universal property that depends only on orientation, 
but also on other intrinsic parameters such as the accretion rate or 
the Eddington ratio. 



5.4 Similarities to non-black hole astrophysical systems 

It is interesting that many apparent parallels can be drawn between 
IRAS 13349, a massive energetic black hole systems fuelling a 
whole galaxy, and other less powerful/smaller systems. For exam- 
ple, the relatively large grains found in the torus of IRAS 13349 
are also encount ered in the dusty winds o f luminous blue vari able 
(LBVs, see e.g. iRobberto & Herbst|[l998l : lo' Hara et al. I I2OO3I) . In 
these extreme massive stars, the maximum grain size is set pro- 
portionally to the mass-loss rate of the high-opacity wind, and 
dust is likely created during stellar eruptions, when this mass- loss 
reaches very high levels (0.001-0.01 AJ,^/vr. lKochanekll20I ll) . In 
IRAS 13349, the high-opacity component (characterised by the 
high optical extinction suffered by the direct light), and the dust 
appear to be co-located, suggesting perhaps also dust formation in 
the IRAS 13349 outflows, be it from the torus or disk. In addition, 
based on sublimation temperature arguments, our placement for the 
iron- grains (detected in the Chandra X-ray spectra) in the rim and 
upper layers of the torus with the silicates (detected in Spitzer IR 
data) in th e interior, gives a pic ture qualitatively similar to that ar- 
rived at bv lChiang et al.l l l200lh for the dust in the surface layers of 
T Tauri and Herbig Ae disks. 
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6 SUMMARY 

Although IRAS 13349 shows many signs of being an obscured 
QSO, the Spitzer spectra show it as a typical QSO. The contin- 
uum peaks at 30 /im; the silicate features at 10 and 18 /im are in 
emission, and PAH features are weak. The emission lines are dom- 
inated by highly ionised species such as [Ne v] and [O iv]. Fur- 
thermore, even though IRAS 13349 exhibits extreme Eigenvector- 1 
characteristics, HET optical spectra argue against its classification 
as a narrow-line-Seyfert 1. The high signal-to-noise and R ~ 1300 
HET data has also allowed us to improve redshift measurements 
for IRAS 13349, locating it at z=0. 10853. 

In the ultraviolet, our low-resolution STIS spectra cover the 
observed wavelength range from 1 150-3 180 A. These spectra show 
for the first time blue-shifted absorption from Lya, N V and 
C IV, with components at systemic velocities of —950 kms~^ 
and —75 kms^^ . The lines are unresolved, and they have intrin- 
sic widths with Doppler parameters that are < 200 kms^^ . As 
seen before (Hines et al. 2001), the UV spectrum is heavily red- 
dened, but at the shortest wavelengths in the UV, below an observed 
wavelength of 1600 A, the spectrum begins to recover in flux and 
appears less heavily reddened. 

Our HETGS spectrum of IRAS 13349 covers X-ray wave- 
lengths from 2-38 A(observed). The continuum is well described 
by a power law with photon index T = 1.9. At low energies we 
observe a soft X-ray excess that can be modelled as a black body 
with kT ~ 100 eV. IRAS 13349 is well known for the ionised 
absorption that is prominent in its X-ray spectrum (Brandt et al. 
1996, 1997; Sako et al. 2001; Longinotti et al. 2003; Holczer et al. 
2007). Our spectrum shows absorption lines attributable to at least 
two warm absorbers at velocities comparable to the higher velocity 
UV absorber. The lowest ionisation WA- 1 is co-located with the 
UV absorber. 

Using simultaneous Chandra HETGS and HST/STIS spectra 
obtained in February 2004, contemporaneous ground-based opti- 
cal spectra from the HET, and archival Spitzer IRS spectra from 
2005 as well as IRAS, ISO, and 2MASS photometry data, we have 
constructed a broad-band view of the SED of IRAS 13349. If we 
compare the observed SED of IRAS 13349 to the median QSO 
spectral energy distribution (SED) of Elvis et al. (1994), we see 
that the IR and X-ray bands are a good match, but the UV and opti- 
cal are severely suppressed. This suggests that these portions of the 
spectrum are either heavily obscured and/or only viewed in scat- 
tered light. By normalising the peak of the infrared continuum and 
the absorption-conected X-ray continuum to the Elvis et al. (1994) 
SED, we determine that the unreddened portion of the optical con- 
tinuum represents a fraction of only 20% of the intrinsic light from 
IRAS 13349. 

If we correct the observed SED of IRAS 13349 for scattering 
and extinction, it peaks at a rest wavelength of ~1000 A. A simple 
sum-of-thermal black body thermal disc spectra provides a good 
fit to the SED. We anive at a black hole mass Mbh ^ 10^ M© 
for IRAS 13349 based on theoretical fits to the UV-Xray SED, and 
the H/3 line width, independently. Theoretical considerations com- 
paring different ionising SEDs reveal that including the UV (i.e. 
disc) as part of the ionising continuum has strong implications for 
the conclusions one would draw about the thermodynamic stabil- 
ity of the warm absorber. Specific to IRAS 13349, we find that an 
Xray-UV ionising SED leads to the conclusion for a continuous 
distribution of ionisation states in e.g. a smooth flow (this paper), 
versus discrete clou ds in pressure equilibrium (previous work, e.g. 
iHolczer et ani2007h . 



To explain the shape of the SED, we developed a geometrical 
model in which we view the nuclear regions of IRAS 13349 along 
a line of sight that passes through the upper atmosphere of an ob- 
scuring torus. This sight-line is largely transparent in X-rays since 
the gas is ionised, but, as previously discussed, it is completely ob- 
scured by dust that blocks a direct view of the UV/optical emis- 
sion region. 20% of the intrinsic UV/optical continuum is scattered 
into our sight line' by the far wall of an obscuring torus. An ad- 
ditional 2.4% of the direct light, which likely dominates the UV 
emission, is Thomson-scattered into our line-of-sight by another 
off-plane component of highly ionized gas. Our model suggests 
that the direct line-of-sight is probably completely obscured by 
dust. In the standard Unified Scheme, such a configuration should 
produce a spectrum dominated by narrow-emission lines, and the 
object would be classified as a Type 2 QSO. However, our model 
also suggests that nearly all of the UV/optical light we see is scat- 
tered from the central engine. Similar inferences have been made 
to explain the observed properti es of o ther infrared-selected QSOs 
— 2MASSI J130005. 3 163 214 jSchmid t et al. 200^ and 2MASX 
J1049433 4-I-5837501 ( Schmidt et al.ll2607l) — and the BALQSOs 
Mrk 231 ( [Gallagher et al.u2005i) . In each case, the view of the cen- 
tral engine, and the scattered light, are significantly obscured and 
the broad-line emission is revealed primarily in scattered light. For 
Mrk 231, imaging polarimetry constrains the scattering region to 
be very compact. For IRAS 13349 we also argue that the scat- 
tering material is fairly close to the BLR. It is likely that viewed 
from a higher inclination, these compact scattering regions would 
also be obscured and the object would be classified as a Type 2 
QSO. In that case, the BLR could be revealed by light scattered 
much farther from the central engine, in an ionisation/scattering 
(bi)cone in direct analogy with other hi ghly polari sed Type 2 QSOs 
(e.g.lHi nes & Wills 1993; Hines et alj 1995. 1999j; lTran et alj200(]|; 
IZakamska et al.ll2005i) . 

Future modeling of IRAS 13349 can take our picture pre- 
sented above as a starting point to calculate more detailed radiative 
transfer models of the scattering processes that include polarisation 
calculations. These can provide more quantitative constraints on 
the viewing angles and sizes of the scattering regions in our com- 
plex view of this intriguing AGN. We note further that few grating- 
resolution X-ray spectra exist for quasar-luminosity objects, and it 
is important to acquire as many as possible so that their proper- 
ties can be compared reliably with those of the well-studied local 
Seyfert galaxies. Similarly, only a few grating-resolution spectra 
exist for strong Fe II, weak [O III] galaxies with soft X-ray spectra 
that are at the extremes of Eigenvector- 1, as is IRAS 13349. 
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